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ABSTRACT 

ON  INVARIANT  SURFACES  AND  BIFURCATION  OF  PERIODIC  SOLUTIONS 
OF  ORDINARY  DIFFEREIMTIAL  EQUATIONS 
by  Robert  John  Sacker 

We  consider  the  ordinary  differential  equation 

dx 

—  =  F(x,u) 

dt 

where  x  is  a  real  n-vector  and  ^  is  a  real  parameter. 
It  is  assumed  that  for  all  n  in  a  neighborhood  of  a 
critical  value  \j.x     this  equation  has  a  periodic  solu- 
tion x=  i{i(t,|a)  which  is  asymptotically  stable  for 
^  <  Hi   (n-1  Floquet  multiples  have  modulus  <1) .   As  n 
increases  through  m  we  assume  a  conjugate  pair  of 
multipliers  \i^)     and  x(^),   leaves  the  unit  disc  caus- 
ing ^     to  become  unstable.   It  is  shown  that  a  two- 
dimensional  asymptotically  stable  invariant  torus  will 
bifurcate  from  ijj  provided  Ca)   x(m)   is  not  a  second, 
third,  or  fourth  root  of  unity,  (b)  d/d^  |A.  (m)  |  >  0,  and 
(c)  a  certain  nonlinear  term  obtained  from  F(x,n)  has 
the  proper  sign  (nonlinear  damping) . 

Some  cases  in  which  (a)  does  not  hold  (resonance 
cases)  are  also  discussed.   In  a  particular  subcase  of  the 
case  \'*(ni)  =  1  (but  x  and  \^  /^  1),  ^     bifurcates 
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into  a  pair  of  subharmonic  solutions,  exactly  one  of 
which  is  stable.  IVhen  a  single  multiplier  leaves, 
xCtii)  =1,   \i;  will  in  general  bifurcate  into  one  stable 
periodic  solution.   The  case  \3(^j_)  =  i   (but  x  /^  1) 
is  discussed  only  by  means  of  an  example  in  which  the 
instability  that  develops  is  so  violent'  compared  to  the 
previous  cases  that  there  is  no  stable  manifold  or 
periodic  solution  in  a  small  neighborhood  of   \i/. 

The  existence  of  the  torus  is  established  by  two 
methods.   The  first  is  to  take  an  n-1  dimensional  hyper- 
plane  S  normal  to  the  orbit   ij;  and  show  that  the  map- 
ping of  S  into  itself,  induced  by  the  vector  field  near 
(ji,  has  an  invariant  curve.   In  the  second  method  the 
torus  is  obtained  as  the  solution  of  a  certain  partial 
differential  equation  which  is  basic  in  the  theory  of 
invariant  surfaces. 

The  second  main  result  is  concerned  with  the  theory 
of  perturbation  of  an  invariant  surface  of  a  periodic 
vector  field.   V/e  consider  the  differential  system 

X  =  f(x,y,(a) 

y  =  gCx^y?!^) 

where  x  and  y  are  real  vectors,   n  is  a  small  real 
parameter  and  the  functions  are  periodic  in  x  (e.g., 
X  represents  the  angular  variables  on  a  torus).  The 
problem  is  to  find  an  invariant  manifold   T(id):  y  =  /(x,^i) 
whenever   t(0)   is  known  and  has  certain  properties.   One 
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of  the  properties  assumed  by  other  authors  is  that  the 
flow  (vector  field)  on  t(0)   is  parallel,  i.e.,  in 
appropriate  coordinates  on  the  surface,  solutions  can  be 
obtained  from  one  another  by  a  translation.   This  restric- 
tion is  dropped  since,  in  general,  parallel  flow  is  lost 
after  a  perturbation  is  made.  The  essential  property  is 
the  relative  strength  of  the  flow  against  the  unperturbed 
surface  and  the  flow  in  the  surface. 

The  function  ;zi(x,n)   satisfies  the  quasilinear 
partial  differential  equation 

)       f  (x,;z^,|a)— —  ~  gCx,/,^)  =  0 

^-~'  ^     ax 

V  V 

which  is  solved  by  an  iterative  scheme.   It  is  linearized 
and  smoothed  so  that  vje  have  a  C*"  equation  Lu  =  h  for 
the  approximate  solution  u.   This  is  made  elliptic  by 
adding  the  Laplacian,   cAu  •*•  Lu  =  h,   c  =  constant,  so 
that  existence  and  smoothness  of  a  solution  becomes  a 
matter  of  quoting  well  known  results.   The  existence  of  a 
solution  of  the  quasilinear  equation  is  obtained  from 
"a-priori"  estimates  for  the  solution  of  the  elliptic 
problem  and  its  derivatives.   If  the  flow  against  the 
surface  is  strong  compared  to  the  flow  in  the  surface, 
then  the  solution  ;e^(x,n)  will  be  very  smooth  since  higher 
derivatives  may  be  estimated. 
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CHAPTER  I 
INTRODUCTION  AND  SUI4MARY 

In  this  dissertation  we  are  concerned  with  the 

ordinary  differential  equation 

dx 

—  =  F(x,n)  Cl.l) 

dt 

where  x  and  F  are  real  n-vectors  and  ^  is  a  real 
parameter.  We  will  say  that   (1.1)  has  an  invariant 
manifold  a,^,   of  dimension  k  <  n,   if  any  solution 
passing  through  a  point  of  c-^     remains  on  a,      for  all 
time  t.  The  manifold  a^  is  said  to  be  asymptotically 
stable  if  all  solutions  of  (l.l)  sufficiently  near  a-^ 
approach  o,   as  t  -  +»,  and  unstable  if  in  any  arbi- 
trarily small  neighborhood  N  of  Oi^  there  is  at  least 
one  solution  which  leaves  N  as  t  -»  +». 

Assume  that  for  |a  <  no  >   a,^  is  asymptotically 
stable,  but  for  ^  >  po   i't  is  unstable.   One  of  the 
simplest  examples  of  this  for  n  =  2,   k  =  0,   is 

2  =  (i+n)z  +  pz^z  (1.2) 

=  i[l+Im  glzj^Jz  +  [n+Re  0|z|^]z 

where  z  =  Xi  +  ixa ,   no  =  0  and  oo  =  0.   The  character- 
istic exponents  for  the  singular  point  z  =  0  are 
a(u)  =  i  +  t-L  and  a.   If  Re  g  <  0  (nonlinear  damping) 

-  1  - 
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then  for  n  >  0  this  equation  has  an  asymptotically  stable 
limit  cycle 


z  = 


\-Re  3/ 

i.e.,  a  stable  invariant  manifold  Oi(u)  with  oi(0)  =  ao » 
which  branches  off  or  bifurcates  from  ao  as  n  increases 
through  zero  (Fig.  1). 

Definition;   W6  say  that  the  mainfold  o   bifur- 
cates*  from  o-^   as  \i     increases  through  no   if 
a   (no)  =  a,-  but  a  C|a)  and  Or.     have  no  points  in  common 
for  n  >  no • 

The  above  example  illustrates  a  general  result 
obtained  by  E.  Hppf  [2]  which  may  be  stated  in  a  somewhat 
weaker  form  as  follows.  Assume  (1,1)  is  analytic  and  has 
a  stationary  solution  x  =  0  for   |[aj   sufficiently  small. 
Assume  that  for  ^  <  0  all  the  characteristic  exponents 
ai[i)      satisfy  Re  a(ia)  <  0.   As  \j.     increases  through  zero 
assume  a  conjugate  pair  of  exponents  aCn)?  aCii)  /  0 
crosses  the  imaginary  axis  and  enters  the  right  half  plane 
with  Re  aCO)  =  0,   Re  a'CO)  >  0?   while  the  remaining 
exponents  satisfy  Re  aCta)  <  s-o      for  some  constant  ao  <  0. 
Then  if  a  certain  constant  depending  on  terms  of  degree  <  3 
in  the  expansion  of  F  about  x  =  0  has  the  proper  sign 
(nonlinear  damping),  then  the  solution  ao*   x  =  0 


See  l"Iinorsky  [1]  for  a  discussion  of  bifurcation, 
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bifurcates  into  an  asymptotically  stable*  periodic  solu- 
tion oi()i):   x=  itiCt,|ji)   for  ^  >  0  sufficiently  small. 
\j((t,0)  =  0  and   \|((t,^)  1^   0  for  la  >  0.   If  the  constant 
has  the  opposite  sign,  then  %     will  exist  in  advance 
(^  <  0)  but  will  be  unstable,  i.e.,  as  \x     decreases 
through  zero  the  unstable  origin  oo  becomes  stable  and 
bifurcates  into  an  unstable  periodic  solution.  Hopf 
assumed  analyticity  of  the  functions  and  used  power  series 
expansions  in  his  proof  —  the  method  of  Poincar^. 

In  this  dissertation  we  consider  the  case  in  which 
ili(t,ii)  becomes  unstable  as  \x.     increases  further  and 
passes  through  a  second  critical  point  m ,   and  a  pair 


of  Floquet  exponents  a(|-i))  aCj-i)   enters  the  right  half 
plane  Re  C  >  0  with  Re  a' Cm)  >  0,   i.e.,  a  pair  of 


Floquet  multipliers  x(^),  \{\x)      leaves  the  unit  disc 
Id  <  1  with  d/d^|\(nl)j  >  0. 

The  Floquet  exponents  are  associated  with  (1.1) 
linearized  along  ^l((t,^Jl), 

y  =  F^(,t(t,n  ),n)y  ^^'^^^i 

If  the  frequency  of     \|ii    s   ijiCtjUi)      is     oui      then 

Im  a(m)   =  uj2      (determined  only  up  to  integral  multiples 


*  i.e.,      n-1     of  the  Floquet   exponents  associated  with 
f(t,|a)     have  negative  real  parts. 
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of  uui)   is  the  frequency  of  oscillation  about  ^x      o^ 

solutions  of  (1.3)    starting  near  iiii .  Floquet  theory 

asserts  that  (1.3)    has  a  family  of  solutions  of  the 

Ml 

form 

j6il)   -   A(t)  cos  u)2t  +  B(t)  sin  ouat 

vhere  A  and  B  have  frequency  uui .   If  cua/iwi   is 
irrational  (non-resonance  case)  then  x^Ct)   is  not 
periodic.   If  lua/wi  =  rational  =  p/q  in  lowest  form 
(resonance  case)  then  /(t)   is  periodic  and  has  frequency 
uui/q,   i.e.,  the  period  P  =  qPi   (Pi  =  2n/u)i   is  the 
period  of  \|(i ) ,   If  q  >  1  v;e  call  /(t)   a  subharmonic 

solution.   In  terms  of  the  Floquet  multiplier   x(|-ii) 

i  2n  ^  a 

=  e     ^^ ,  non -resonance  is  characterized  by  X  (m)  /  1 

for  any  integer  q  and  resonance  by  X  (i-ii)  =  1  for  some 

integer  q. 

Throughout  this  work  we  will  assume  that  F  has 
p  >  5"  continuous  derivatives  with  respect  to  x  and  \j, 
in  some  neighborhood  of   IjI(t,^).   V/e  will  prove  the  follow- 
ing. 

Statement  1;   If  certain  low  order  resonances  are 
avoided,  specifically,  if 

(a)  uua/wi  ^   p/q  in  lowest  form  with  q  =  1,2,3,1+ 
i.e.,  if  X*^(m)  /  1  for  those  q,  and  if 

(b)  Re  a'(>^i)  >  0  and 

(c)  a  certain  constant  obtained  from  the  terms  of 
degree  <  3  in  the  expansion  of  F  about   i)((t,ni)  has 
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tr.e  proper  sign,  then  yCtju)  birurcates  into  a  unique 

asymptotically  stable  torus  aaC^)  as  \i   increases  through 

m.  For  r  an  integer,   0  <  r  <  p-1,  02  (n)   exists  as  an 

r -times  continuously  d liferent iable  manifold  for 

0  <  ^a  -  (ij.  <  6   where  6   depends  on  r,   the  terms  of 

degree  <  3  mentioned  in  (g) ,  and  Mo  =  [Fj  *.   In  general 

6  -♦  0  as  r  -  00.   If  the  constant  in  (c)  has  the  opposite 
r 

sign,  then  the  torus  02     exists  in  advance   (la  <  ni)  and 
is  unstable. 

A  model  of  such  a  bifurcation  phenomenon  was  given 
by  E.  Hopf  [3],  Hence  we  show  essentially  that  for  the 
second  bifurcation  (i.e.,  periodic  solution  into  torus) 
Hopf's  model  is  the  generic  situation. 

It  is  worth  noting  that  even  if  F(x,u)   is  analytic, 
this  torus  may  not  be  analytic  but  only  represented  by 
functions  with  finitely  many  derivatives.  Thus  any  method 
based  on  power  series  expansions  must  fail. 

Three  resonance  cases  are  discussed  and  we  prove 
the  following. 

Statement  2;   If 

(a)  tua/oui    =  IA5      i.e.,      x'^Cm)    -  1 

(b)  Re   a'    Ciai)    >  0  and 


*  This  norm  is  defined  at  the  end  of  the  chapter.  The 
differentiation  is  with  respect  to  x  and  |a  and  the 
supremum  is  taken  over  a  small  neighborhood  of  \(iCt,ia). 
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(c)  certain  conditions  are  satisfied  by  the  terms  of 
degree  <  3  in  the  expansion  of  F  about  ^jiCtjm),   then  as 
n  increases  through  m,  ili(t,fa)   bifurcates  into  a  pair  of 
subharraonic  solutions.*  One  solution  is  asymptotically 
stable  and  the  other  is  unstable. 

The  next  result  deals  with  the  case  in  which  one  real 
Floquet  exponent  ai\s)     becomes  positive. 

Statement  3:   If 

(a)  a(0)  =  0,   a(>^)   real,     (b)  a'C^)  >  0    and 

(c)  a  certain  constant  obtained  from  the  terms  of  degree 
<  2  in  the  expansion  of  F  about  ij)(t,m)is  not  zero,  then  as 
H  increases  through  ui  ,  i|((t,|i)  bifurcates  into  an  asymp- 
totically stable  periodic  solution  having  the  same  period.* 

The  third  resonance  case,  wa/uui  =  1/3 »   is  discussed 
by  means  of  an  example  which  illustrates  the  violent 
instability  which  develops  after  bifurcation. 

The  resonances  (Statements  2  and  3  and  the  example) 
are  treated  in  Chapter  III  where  the  rate  of  growth  of  the 
bifurcating  manifolds  is  also  discussed.  The  proofs  of 
Statements  2  and  3  require  the  usual  construction  of 
periodic  solutions  by  means  of  the  implicit  function 
theorem  Ih],      In  contrast,  Statement  1  requires  the  con- 
struction of  a  2-dimensional  surface  and  thus  different 
techniques  are  needed.   Statement  1  is  proved  in  Chapter  II 
and  again  in  Chapter  III.   In  Chapter  II  the  problem 


*  The  deviation  of  period  with  n  is  eliminated  by  intro- 
ducing a  new  parameter  s  in  place  of  t. 
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is  treated  by  working  vfith  the  mapping  induced  by  the 
vector  field  near   iji,   the  so-called  "method  of  surfaces 
of  section"  of  Poincar^  and  G.  D,  Birkhoff .  The  method 
is  described  in  detail  in  the  introduction  to  that  chap- 
ter. 

In  Chapter  III  we  treat  the  differential  equation 
directly.  The  torus  is  obtained  as  the  solution  of  a 
certain  partial  differential  equation  which  arises  in  the 
theory  of  perturbation  of  an  invariant  surface  of  a 
periodic  vector  field.   This  theory,  treated  in  Chapter  IV, 
is  concerned  with  the  problem  of  finding  an  invariant  mani- 
fold t(u)  of 

X  =  f  (x,y,ia) 

Cl.^) 
y  =  g(x,y,^) 

whenever  t(0)   is  known  and  has  certain  properties,   n 
is  a  small  parameter,  x  and  y  are  vectors,  and  the 
functions  are  periodic  in  x.   If  t(0)   is  the  surface 
/o (x) ,  then  the  problem  is  to  find  a  solution  tCh); 
y  =  /(x,n)   of  the  quasilinear  partial  differential  equa- 
tion 

f  (x,/,;a)  — ~  -  gU,j6,\A)    =    0        (1.5) 
^       dx 

V  V 

with  /$(x,0)  =  /oCx).  The  properties  of  the  flow  (vector 
field)  in  the  surface  tCO)   and  the  flow  against*  the 


*  We  do  not  consider  the  case  in  which  the  flow  is  away 
from  t(0). 
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surface  are  essential  in  determining  the  existence  and 
smoothness  of  the  perturbed  surface.  We  say  that  the 
flow  in  t(0)   is  parallel  if  x  =  constant  on  the  surface. 
The  flow  is  parallelizable  if  it  is  parallel  in  appropriate 
coordinates  on  the  surface.  The  case  of  parallel  flow  has 
received  much  attention  (Diliberto  [^,  6],  Bogoliubov  and 
Mitropolski  [7,  8],  Hale  [9,  10],  Diliberto  and  Huff ord 

[11]). 

In  treating  the  general  problem  of  perturbation  of 

invariant  surfaces  the  assumption  of  parallel  flow  is  much 
too  stringent  and  unnatural  since  it  is  generally  lost 
after  the  perturbation  is  made.  When  the  unperturbed  sur- 
face contains  lower  dimensional  invariant  manifolds  (e.g., 
stationary  points)  it  is  usually  impossible  to  parallelize 
the  flow.  The  characteristics  on  the  surface  may  form 
envelopes  and  it  seems  unsatisfactory  to  use  methods  such 
as  those  of  Bogoliubov  and  Mitropolski  [7]  in  which  (1.5) 
is  solved  by  integrating  along  characteristics.   The  char- 
acteristics method  is  associated  with  initial  value  prob- 
lems whereas  the  problem  of  solving  (1,5^)  defined  on  a 
torus  may  be  formulated  more  naturally  as  a  boundary  value 
problem  in  which  the  boundary  condition  to  be  imposed  is 
the  periodicity  of  the  solution  in  x.   It  is  the  latter 
formulation  which  is  considered  in  Chapter  IV.   The  solu- 
tion of  (1.5")  is  obtained  by  an  iterative  process.   At 
each  stage  the  equation  is  linearized  and  smoothed  so  that 
we  obtain  an  equation  Lu  =  h  where  h  and  the 
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coefficients  in  L  are  c"  functions.  This  is  then 
modified  to  give  an  elliptic  equation 

cAu  +  Lu  =  h 

where  c  =  c(n)   is  a  constant  which  tends  to  zero  as  n, 
the  number  of  iterations,  tends  to  oo  and  a  is  the 
Laplacian.  Under  the  given  conditions  we  obtain  existence 
and  C"  smoothness  of  a  solution  of  the  above  equation. 
The  existence  of  the  solution  of  the  nonlinear  equation 
then  follows  from  the  "a  priori"  estimates,  fbr  the  solution 
of  the  linear  equation  and  its  derivatives.   In  deriving 
these  estimates  there  are  two  parameters,   g  >  0  and 
X  <  0,  which  play  a  fundamental  role,   p  is  a  measure  of 
the  flow  against  the  unperturbed  surface  and  x     is  a 
measure  of  the  flow  in  the  surface.   It  is  shown  that  if 
3  +  rx  >  0,  r  an  integer  >  2,  and  another  condition*  is 
satisfied,  then  for  sufficiently  small  n  the  perturbed 
surface  t(h)  has  r-l  Lipschita  continuous  derivatives. 
This  is  the  second  main  result  established  in  this  disser- 
tation (Theorem  8),   It  says  essentially  that  if  the  flow 
against  the  surface  is  great  (  3  large)  compared  to  the 
flow  in  the  surface,  then  the  perturbed  surface  will  be 
very  smooth.  Heuristically  we  can  see  the  need  for  such 


*  This  condition  (number  3  in  Theorem  8)  is  used  to  solve 
the  linearized  equation  and  is  precisely  the  "positive" 
condition  of  Friedrichs  [12]. 
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an  assumption.  For  consider  (in  one  x  dimension)  two 
characteristics  C^   and  Gg  approaching  the  equilibrium 
X  =  0,  (see  Figs.  2  and  3  where  t(0)   is  the  line  y  e  0.) 
If  as  in  Fig,  2,  the  flow  against  y  =  0  is  small  com- 
pared to  the  flow  in  y  =  0  (shown  by  arrows),  then  a 
cusp  may  develop.   McCarthy  [13]  and  Kyner  [1^]  discuss 
an  example  of  such  a  situation.   If  the  relative  strength 
of  these  flows  is  reversed  (Fig.  3)  then  C^  and  Cj 
tend  to  meet  in  a  smooth  fashion. 

Thus  even  if  f  and  g  are  C°°  functions  and 
t(o)   is  a  C"  surface  it  may  happen  that   T(|a)  has  only 
finitely  many  derivatives.   In  fact,  in  the  example  of 
Kyner  the  perturbed  surface  has  a  cusp  for  arbitrarily 
small  n  >  0. 

The  result  of  Kyner  [IM]  is  not  restricted  to 
parallel  flow.   In  the  case  of  non-parallelizable  flow 
he  postulates  the  invertlbility  of  a  certain  linear 
transformation  and  obtains  a  perturbed  surface  which  is 
Lipschitz  continuous.  Higher  derivatives  are  not  treated. 

In  the  special  case  in  which  the  flow  on  t(0)   is 
parallel,   X  =  0  and  thus  rda)  will  have  any  finite 
number  of  continuous  derivatives  provided  n  is  suf- 
ficiently small.   Such  is  the  case  if  we  consider  the 
equation 

i  =  F(§)  +  nFn5,t,u)  (1.6) 

(§  an  n-vector)  where  the  functions  are  C°°   in  all 
arguments  and  have  period  T  in  t.   It  is  assumed  that 
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Fig-.  I 


rco) 


Fi  G-.    3 


rc) 
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for  n  =  0  there  is  a  limit  cycle  r  which  is  asymp- 
totically stable  (n-1  of  the  characteristic  exponents  of 
the  linear  variational  equation  have  negative  real  parts). 
tCO)   Is  then  the  torus  rX{0  <  t  <  Tj   with  the  faces 
t  =  0  and  t  =  T  identified  (see  Fig.  h) .      It  is  known 
([1^],  [8],  [9],  [11])  that  coordinates  can  be  intro- 
duced such  that  the  flow  on  t(0)   is  parallel.  Hale  [9] 
and  Dilibetto  and  Hufford  [11]  treat  (1.6)  in  the  degen- 
erate case  in  which  the  flow  against   t(0)  vanishes  for 
H  =  0  (some  of  the  characteristic  exponents  have  zero 
real  part  for  n  =  0) .  It  is  this  case  which  arises  in 
the  bifurcation  problem  of  Chapter  III  (see  equation  3 '2^) 
and  thus  the  results  of  these  authors  could  have  been 
applied  at  this  stage  to  prove  existence  of  the  torus. 
The  theory  of  Chapter  IV  also  takes  this  degeneracy  into 
account  and  may  be  considered  as  a  new  derivation  of  the 
results  of  the  above  authors.   Its  chief  advantage,  how- 
ever, lies  in  its  applications  to  cases  of  non-paralleliz- 
able  flow.   (see  the  example  in  Chapter  IV). 

Notation;   V/e  nov/  introduce  some  notation  which 
will  be  used  throughout  the  rest  of  our  work.  For  r  a 
positive  integer 

F(x,y,n)  €  C^(x,y)  n  C(^) 

means  that  F  has  r  derivatives  with  respect  to  x 
and  y  which  are  continuous  in  (x,y,ia). 
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In  this  figure  r  is  assumed  to  be  in  the   (§i»l2)  plane. 
The  coordinates   (§3 »•••»§  )  a^©  1^°^  shown. 
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F(x,y,^)  6  C^(x,y)  n  Lip^(x,y)  n  C(|a) 

means  the  same  as  above  and  the  derivatives  are  Lipschitz 
continuous  in  x  and  y  uniformly  in  la. 
For  a  vector  f  =  (f.)  define 

If  I  =  (Z  f?)   =  (f,f)^ 
i  ^ 

We  will  consider  functions  f Cx,n)  defined  for  x  a 

vector  in  some  domain  G  and  |a  a  small  parameter.  For 

r  >  0  an  integer  define  the  norm 

l^lr  =  0|?|r  l^g  P'f| 

where  D^  is  any  p    order  derivative  of  f  with 

respect  to  the  components  x.   of  x.  Unless  otherwise 

stated J  differentiation  will  never  be  performed  with 

respect  to  \i. 

For  two  real  vectors  u  =  Cui,,.,,u^)  and 

V  =  (vi,,.,,v  ),   the  symbol   (u,v)  will  denote  the 

usual  inner  product 

n 

(u,v)  =  >   Uj^V^   . 


CHAPTER  II 
BIFURCATION  -  MAPPING  METHOD 

1.  Introduction 

In  this  chapter  \ie   treat  the  bifurcation  problem 
by  considering  the  mapping  induced  by  the  vector  field 
near  the  periodic  solution:   the  method  of  surfaces  of 
section  mentioned  earlier. 

We  consider  the  differential  equation 

dx 

—  =  FCxju)  (2.1) 

dt 

where     x     and     F     are  real  n-vectors  and     ^     is  a  real 
parameter.      Suppose  that   for  all  sufficiently  small     p, 
llj]   <  ^   ,      (2.1)   has  a  periodic   solution     x  =  ijiCtjia)   with 
period     2n.      Assume  that     F     has     p   >   5"     continuous 
derivatives  with  respect   to     x     and     n     in  some  neighbor- 
hood of     \|)     and  that   in  this  neighborhood 

tFlp<Mo 

where  differentiation  is  with  respect  to  x  and  \j.. 
Suppose  that  ^     is  asymptotically  stable  for 
H  <  0  in  the  sense  that  n-1  of  the  Floquet  multipliers 
lie  inside  the  unit  circle  in  the  complex  plane.   As  n 
increases  through  zero  we  assume  that  a  con;)ugate  pair 
of  the  multipliers 
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2TTa(!-i)        2TTa(|i) 

\C|a)    =  e  ,      x(n)    =  e 

leaves  the  unit   circle  with 

Re   a(0)    =  0,     Re  a'(0)   >  0 

while  the  remaining  n-3  stay  inside.  Hence  the  periodic 
solution  iji  becomes  unstable. 

We  do  not  require  that  ij     is  obtained  by  a  bifur- 
cation from  an  equilibrium  as  discussed  in  Chapter  I  but 
just  consider  any  periodic  solution  which  loses  its 
stability  in  the  above  manner  as  [i     passes  through  a 
critical  value  (assumed  to  be  n  =  0). 

2.  Reduction  to  a  Mapping  Problem 

We  now  obtain  the  mapping  described  by  the  flow 
near  the  periodic  solution  ^.     Let  C  be  the  curve 
described  by  ^.      Choose  a  point  P  on  C  and  let  S  be 
the  n-1  dimensional  hyperplane  normal  to  C  at  P.   Let 
P  be  the  origin  of  coordinates  x  =  ^^(  1") » •  •  •  >-^(n-l)^ 
describing  S.   It  is  known  [16;  p.  ^0]  that  any  solution 
of  (2,1)  starting  on  S  and  close  enough  to  P  will 
return  to  S  for  the  first  time  after  a  time  lapse  of 
2tt  +  e  (x)   where   e  -  0  as  x  -  0.   In  a  small  neighbor- 
hood of  P  this  defines  a  mapping  T  of  S  into  itself 
having  P  as  a  fixed  point 

T:   Xi  =  Dx  +  ....  (2.2) 
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It  is  kno-wn  [16;  p.  I63]   that  the  eigenvalues  of  D  are 
just  the  n-1  previously  mentioned  Floquet  multipliers. 
Thus  for  n  <  0,   P  is  an  asymptotically  stable  fixed 
point  which  becomes  unstable  as  [i     increases  through  zero. 

The  existence  of  a  torus  will  be  established  by 
showing  that  P  bifurcates  into  a  closed  curve  C,  lying 
in  S  which  is  invariant  under  T^i.e.,   if  q  is  a 
point  of  C,   then  Tq  =  qi   is  on  C   (see  Fig.  ^) . 


3.  Normal  Form  for  the  Mapping 

By  appropriately  transforming  coordinates  we  will 
put  the  mapping  in  a  certain  normal  form  which  is  needed 
in  the  existence  theorem.   We  assume  that  a  real  linear 
transformation  of  coordinates  has  been  carried  out  so  that 
D  is  in  the  real  form 


D  = 


Re  \(n) 
Im  \(|a) 

-Im  \(^)   ' 
Re  \(n) 

0 

0 

S(^) 

where  S(n)   is  a  matrix  whose  eigenvalues  o(n)   satisfy 

la(0) I  <  1  (2.3) 


and 


sup   |SCO)Y|  =  ao  <  1 
1Y|=1 

where  Y=  Col  ^x^^,-) » . . .  ,X/  ,  s )  and   |Y|   is  the  ordinary 
Euclidean  length  of  a  vector   (YjY)"^, 
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Fig-.  5" 
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Let     z  =  X-    ,+  ix.    > .     Then  in  the        z,Y     coordinates  the 
mapping  is 

Zi   =  x(n)z  +  U(z,z,Y,^) 
T:  iZ.h) 

Yi    =  SCn)Y  +  VCz,z,Y,n) 

vhere  U  is  complex  and  V  is  real.  The  functions 
U,  V,  X  and  S  have  five  continuous  derivatives  in  a 
neighborhood  of   (z,Y,|i)  =  0  which  we  assume  to  be 


zr  "    1^1  ^  2,     |nl   <  u 


« 


We  assume  U  and  V  to  be  expanded  into  Taylor  series 
up  to  polynomials  of  degree  h   plus  remainder. 

In  order  to  motivate  the  following  theorem  consider 
the  following  example  of  a  mapping 


2 


i+^+3|z| 
Zi  =  e        z 


Yi  =  SY  H-  b]  z| 


2 


where  3  is  complex,  b  is  a  real  constant  vector  and 
S  is  a  matrix  which  satisfies  (2.3) •   If  Re  g  <  0,   this 
mapping  has  the  invariant  curve 

/      i 

,   Y  =  CS-I)   b 


\-Re  g/  .  -Re  p 

contained  in  a  neighborhood   |z|  <  c^M,   |Y|  <  c^n  ,   c 
a  constant.  We  now  consider  the  effect  of  perturbing 
this  example  by  adding  more  nonlinear  terms.   In  such  an 
oblate  neighborhood  Y  is  small  of  order  n  whereas  z 
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is   small  only  to  order     \i   .      More  precisely  consider  the 
neighborhood 

N:      z  =  aC,        1=3^^ 

where   |Cl^  +  \'^f  ^   2  and  a  Is  a  small  parameter, 
0  <  a  <  1.   In  N  the  monomial  M  =  Y?.  >,  2*^2  , 
T  =  2p  +  q  +  r,   satisfies   |M  L  <•  c(k)a''^  where  differ- 
entiation is  with  respect  to  ^(j_)>  C  an<l  C  an<i  c^^) 
is  a  constant  depending  only  on  k. 

Definition.   We  call  t  the  weight  of  M^.   In 
the  following  theorem  we  will  transform  the  mapping  (2A) 
into  a  form  similar  to  the  above  example  in  which  monomi- 
als of  weight  2  and  3  will  serve  to  determine  the 
invariant  curve  in  the  first  approximation  while  the 
remaining  terms  will  act  as  small  perturbations.  The 
transformations  used  are  of  the  type  treated  by  B,  Segre 
[17].  See  also  A,  Kelly  [18]  for  an  excellent  bibliog- 
raphy. 

Theorem  1.   If  x^(0)  /  1,  \3(o)  /  l  and  (2.3)  is 
satisfied,  then  there  exists  a  no  <  M«  such  that  for 
\\j.\      <  uo   there  exists  a  transformation 

z  =  cw  +  P(w,w,W,c,^i) 

0  <  c  <  1 
Y  =  cV/  +  P(w,w,W,c,ti) 

,  with  P  and  P  polynomials  in  w,w,  and  W,  which 
carries  C2.M-)  into  the  form 
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2TTaCu)   +  c23(^)lwl2 

Wi    =  e  w  +  Ru.(v;,w,W,c,n) 

(2.5) 
Wi    =  S(n)W     +     R3(w,vJ,W,c,n) 

where     R      are   function  whose  Taylor  expansions*  about 
(wjWjW)    =  0     contain  no  terms   of  weight   <  k,     P     and     P 
contain  no  constant  or  linear  terms.     For  a  proper  choice 
of     c     the  mapping  is   one-to-one   in  the  neighborhood 
|w|2  H-    |W|2   <  2. 

Proof;      In  this   proof   subscripts   on  functions  will 
carry  the  same  significance  as  for     R,      defined  above.      We 
write   the   first   equation  of    (2,^-)    as 

zi    =  xCu)z  +   u(z,z,ia)    +  v(z,z,Y,|i)    +  Fi^  (2.6). 

where  u  contains  quadratic  and  third  order  terms  in  z 
and  z  and  v  =  Y'[p(ia)z  +  q(|a)z]  with  p,  q  vectors 
and  prime  denoting  transpose.     We  transform   (2.6)   by 

z  =  w  +  r(^)w^w'^,  k  +  -t  =  2,   3     (2.7) 

to  obtain 

-  k~^  ~  ~ 

wi    =  x(|a)w  +  u(w,w,(i)    -  g(iJi)w  w     +  v(w,w,Y,^)    +  Fi^ 

where 

g(n)    =    Cx^X"^-  X)r(|a),  \   =  \(|a). 

Unless  k  =  2,  ^  =  1  we  see  that  rCia)  may  be  determined 
so  that  (2.7)  removes  the  term  w  w'^  from  u.  By  succes- 
sive applications    of   (2.7))   we   obtain  the  form 


By  the  smsothness  assumption  made  after    (2.V),     R,      is 
well  defined  provided     k  <  ^. 
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wi    =  Xi^)\i  +  b(|a)|w|2w  +  vCw,w,Y,|a)    +  Qu.        (2.8) 

The  transformation 

w  =   C  ■*■  Y'rC^L)C^C^,  k  +  ^  =  1  (2.9) 

with  ^(^)  a  vector  carries  (2.8)  into 

Ci  =  \(^a)C  +  b(n)|cl^C  +  v(c,C,Y,li)  -  Y'g(n)C^C^  +  Gv 

where 

g(n)  =  cs'(^)x^x'^  -  xi]r(u),   X  =  x(^i) 

and  I  is  the  identity  matrix.  Since  X.(0)   is  on  the 
unit  circle  and  from  (2.3)  the  eigenvalues  of  S'(0)  are 
inside  the  unit  circle,  we  see  that  r(|i)  may  be  deter- 
mined for  ^  sufficiently  small  such  that  the  term 

k  t 
C  C   is  removed  from  v.  The  mapping  now  has  the  form 

Ci  =  \(n)C  +  b(^)lc|2c  +  Hi,        (2.10) 
The  second  equation  of  (2."+)  may  be  written 

Yi  =  S(ia)Y  +  sCc,C,n)  +  H3  .         (2.11) 
where  s  is  quadratic  in  Q     and  Q.     The  transformation 

Y  =  X  +  r(n)c^C^,     k  +  ^  =  2 
reduces  (2.11)  to 

Xi  =  S(|a)X  +  s(c,C,n)  -  g(M)C^C'^  +  H3 


where 


g(u)  =  [\^\^I-S(ia)]r(M),      X  =  X(u) 


Just  as  before  r(|a)  may  be  determined  for  la  sufficiently 
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small.  Thus  the  mapping  has  the  form  (2.10)  together 
with 

Xi  =  sc^)x  +  V3 
Using  xCm)  =9  2^aC^i)^   (2.10)  may  be  written 

^,  =  e^'^'^^^hl   -^  3(n)|C|^]C  +H^ 
=  32na(M)+p(^)lCr^  ^H^ 

with  g  =  b/x .   All  the  previous  transformations  may  be 
combined  into  one;   z  =  C"^«'»  jY  =  X+..,  where  the 
dots  represent  polynomials  in  CjCjX  without  constant  and 
linear  terms.   It  is  on-to-one  in  a  neighborhood 
|Cl^  +  jX|^  <  2c.  The  transformation  c  =  cw,  X  =  cW 
then  gives  the  desired  results. 

h.   Existence  of  an  Invariant  Curve 

Having  the  mapping  in  the  form  (2.5)  we  may  now 
establish  the  existence  of  an  invariant  curve  and  hence, 
according  to  the  discussion  of  Section  2,  an  invariant 
torus  of  the  differential  equation  (2.1). 

Theorem  2.   In  (2.H-)  assume 

1.  X^(0)  /  1,  x^(0)  /^   1   (non-resonance) 

2.  x(n),  S(m),  U,  V6C'^(z,z,Y,u)   for  I     an 
integer  >  5,  ||j|  <  n^^,  jzj^  +  |Y]2  <  2   and  assume  S(^) 
satisfies  (2.3).. 

In  (2,5)  assume 

3.  A  =  2n  Re  a'(0)  >  0,   B  =  Re  p(0)  <  0. 
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Let  r,   l<r<'C'  be  a  positive  integer.  Then 
there  exists  a  <  [i»     such  that  for  0  <  n  <  ij_)   (2.^) 
has  an  asymptotically  stable  invariant  curve 

z  =  aoTiIe   +  |ar(e,(j)e 

where  ao  =  J   -A/B.   f  and  g  are  defined  for  all  Q 
and  0  <  ^  <  n  ,  have  period  2n  in  9, 

f,g€C^"Ne)  n  Lip^~^(e)  n  c(n) 

and 

Iflp.i,  |g|r.,  <  c(r) 

where  2(r)   is  a  constant  which  depends  only  on  r  and 
differentiation  is  with  respect  to  9  only.   In  general 
n  -  0  as  r  increases.   ^   depends  on  the  coefficients 
of  terms  of  degree  <  3  in  (2.H),  i.e.,  on  terms  of  degree 
<  3  in  the  expansion  of  (2.1)  about   iii.   u   also  depends 
on  the  constant   M©   defined  after  (2.1). 

Proof.   In  (2.5^)  we  assume  c  to  be  fixed  such 
that  the  mapping  is  one-to-one  for   ]wj^  +  |W|^  <  2. 
Note  that   c  depends  only  on  the  coefficients  of  terms  of 
degree  <  3  in  (2. ^) . 

We  restrict  attention  to  a  neighborhood  w  =  az, 
W  =  a^Y  (2.13) 

where  \z\^   +  [Yj^  <  2  and   0  <  a  <  1.   Then  (2.5')  becomes 

2TTa(n)+a2c2p(n)  |zl2 
Zi  =  e  .  z  +  fi (z,z,Y,a,^) 

Yi  =  S(^)Y  +  f2(z,z,Y,a,ia) 
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and 


fiL  <  c(r)a3,     if2!   <  cCr)a 


where  differentiation  is  with  respect  to  z,z  and  Y. 
Note  that  the  subscripts  on  the  functions  are  used  simply 
for  identification  and  have  nothing  to  do  with  "weight" 
as  they  did  in  the  previous  theorem. 

Expanding  a,  3,  and  S,   ye  obtain 

Zi  =  exp[2TTaC0)  +  ^2TTa'(0)  +  a^c^gCO)  jz|  ^Jz  +  f3  (z,z,Y,a,n) 

Yi  =  S(0)Y  +  f;^(z,i,Y,.a,n) 

where 

jfgl  <  c(r)(M2+a2n+a3),    jfi,]  <  c(r)(a+u). 


Letting 


S(0)  =  S,       2na(0)  =  ij 
2na'C0)  =  A  +  iA,     p(0)  =  B  +  iB 


and  choosing 


An 


c^B 


c^B 


0  <  la  <  ^10    =  minCia    ,  ) 

*     -A 


we  obtain 


Zi    =  expCi(T+p/)    +  nA(l-|z|2)]z  +  f5(z,z,Y,^)      (2.1^) 
Yi    =  SY  +   f6(z,i,Y,|a) 


where 


and 


AB 

/  =  A |z|2 

B 


3 


jreL  ^  c(r)n 


J. 

2 
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The  unperturbed  mapping   (fj,!^  =  0)  has  an 
invariant  curve   jz]  =  1,  Y  =  0.  We  introduce  co- 
ordinates in  a  small  toriod  (see  Fig.  6)  surrounding 
this  curve. 

■5fl  1 

z  =  (1  +  p)e^'',   1p1^  +  lYj^  <  -  .       (2.1?) 

The  first  equation  of   (2.m-)   becomes 

Pi    =   (l-2uA)p    -  3iaAp2    _  |_^^p3   +  f7(p,e,Y,^) 

AB 

©1    =  9  +   T   +  ia[A (1+p)^]   +  f8(p,e,Y,n) 

B 

where 

3 

Letting     p=iau,        Y  =  |av,     v;e  finally  obtain  (2.16) 

Ui    =(l-2nA)u  +     ki[g(©)    +  G(e,u,v,^)] 

(2.17) 

©1     =    e    +    T  +     u[fo      +    F(©,U,V,u)] 

Vi    =  Sv  +     h(e)   +  H(G,u,v,^) 

where     fo      is  a   constant,   all  functions  have  period  2tt 

in     ©     and 

lEl^,    |h|^<c(r);    |GI^,    1H|^,    1F|^  <  c(r,K)^        (2.18) 

for 

1 
juj^   +    ivl^  ^  kS        0  <  h  ^  

II  I       I        -»  -    r-    --    ^^2 

where  differentiation  is  with  respect  to  all  variables 
except  u.   c  and  c  are  constants  which  depend  only  on 
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Ft  0"   G 


1  - 


=  SCAIAK 


In  each  figure  the  toroidal  neighborhood  is  obtained  by- 
rotating  the  circle  about  the  vertical  axis. 


ri  0-'  7 


Vv/   =  SC^LAR 


\^^l 
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the  arguments  shown.  See  Fig.  7  for  a  description  of  the 
toroidal  neighborhood  in  the  W,v;  coordinates.  Finding 
an  invariant  curve  u  =  u(9,|a),  v  =  v(9,|j)  of  (2.17)  is 
equivalent  to  solving  the  functional  equation 


(2.19) 


uOi)  -  (l-2nA)u(e)  =  |aG(0,u,v,u) 

v(9i)  -  sv(e)     =  H(e,u,v,n) 

01  =  e  +  T  +  |aF(9,u,v,ii) 

where  F  =  fo  +  F,   fo  a  constant,   G  =  g(9)  +  G  and 
H  =  h(9)  +  H  satisfy  (2.18).  Theorem' 3  which  follows 
guarantees  the  existence  of  a  unique  solution  u  =  u(9,|i), 
V  =  v(9,[a)   of  (2.19)  satisfying  all  our  requirements. 
The  form  (2.12)  is  obtained  by  applying  the  transformations 
(2.16),  (2.15),  (2.13)  and  finally,  the  transformation  of 

Theorem  1. 

Vie  now  prove  the   stability  of  the  invariant   curve. 

If     u(9,^),    v(9,|i)      is   the   solution  of    (2.19)    for 

0  <  la  <  (I  let     u  =  u(e,(i)    •+  6'u,    v  =  v(9,|i)    +  6v     in 

(2.17)  to  obtain  the  nonlinear  variational  mapping 

6U1  =  (l-2nA)6u  +  u(G^~u'F^)6u  +  m(G-^-u 'F^)  6v 

6Vi  =  S6V  +  (H^-uv'F^)6u  -»■  (H^-^v'F^)6v 

where  prime  is  d/d9  and  F^  etc.  are  evaluated  at 
u(9,^)  +  efiu,  v(9,|i)  +  e6v,   0  <  e  <  1.   From  (2,l8)  we 
see  that  the  derivative  of  F,  G,  and  H  with  respect 
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to  u  and  v  are  <  (constant)  n"^  uniformly  in  a  small 
tube  R:   |6u|2  +  ]5v|^  <  constant.  Hence  in  R  for  ^ 
sufficiently  small,  0  <  n  <  |j  <  n  ,  the  above  mapping  is 
a  contraction.  This  completes  the  proof  oi"  the  theorem. 

^.  Solution  of  Functional  Equation 

The  above  stability  argument  is  based  on  the  prin- 
ciple of  contraction,  i.e.,  the  eigenvalues  of  the  linear 
part  of  the  mapping  are  less  than  unity  in  modulus.  This 
same  principle  will  now  be  utilized  to  prove  existence  of 
a  solution  of  equation  (2.19). 

Theorem  ^.   In  (2.19)  assume 

1.  A  >  0,   |Sv|  <  ao|v|,   0  <  ao  <  Ij   oo   constant 

2.  F,  G,  H  6  C^(e,u,v)  n  C(iji),   r  >  1, 

for  0  <  ^  <  ^10  and  |u  |  ^  +  [v  j  ^  <  K^ ,  K  sufficiently 
large  and  assume  these  functions  satisfy  (2,l8)  and  have 
period  2tt  in  9. 

Then  there  exists  n  <  ^q      such  that  for 
0  <  n  <  M   (2.19)  has  a  unique  solution  u(9,n),  v(0,n) 
having  period  2Tr  in  0. 

u,v  €  C^'^iQ)    n  Lip^"^(9)  n  C(n) 

and   lulj,_j^,  lv|    <  K  uniformly  in  |j  where  differ- 
entiation is  with  respect  to  0  only,   np  depends  on 
A,  00,  r  and  the  constants  c  and  c  in  (2.18).   In 
general,  n^  -  0  as  r  increases. 

Proof.  Define  b  =  min(2A,l-ao ) .  Then     0  <  b  <  1. 
Also  define 

K  =  b    Co  Cr)  c(r)  2 
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where  Co(r)   and  6  are  defined  in  Lemma  2  and  c(r) 
comes  from  (2.18).  We  will  construct  successive  approxi- 
mations u   and  V   with  uo  =  0,  vo  =0.  For  con- 
n       n            ' 

venience  let  w  =  CiJ^).   As  an  induction  assumption  suppose 

n    Vj^ 

w     satisfies 
n-i 

■"n-i  ^  ^^^®^  ^  ^^^^'     0  <  n  <  uo      (2.20) 
l^n-Jr<^  ^2.21) 

Inserting  this  approximation  in  (2,19)  we  obtain 


uOi)  -  (l-2nA)u(e)  =  nG^(©,u) 


vOi)  -  sv(e)     =  H^(e,u) 


(2.22) 


01   =  9  +  T  +  iaPn(e,u) 
to  be   solved  for     w   (9,fi),     where 

F^(e,n)  =  F(o,w^^^(e,n),  ^)       (2.23) 

and  so  on.   We  now  verify  the  conditions  of  Lemma  1.  * 

We  have  F^,  G^,  H^  e   C^O)  n  C(^)  for  0  <  ^  <  n'o 

=  min(^0)  IAk^).   Also  if  we  define  m  =  inf^   F^   we 

^  W?U      hQ 

have,  using  (2.I8),  b  +  rm  ^  b  -  cr^^  >  1^/2     for  ^ 
small  ,0<n<^'<ia'o,  where  c  depends  on  K  and 
c.(r,K).  Unless  F  =  0  we  see  that  n'  -  0  as  r 
increases.   By  Lemma  1,   (2.22)  has  a  solution  w  (0,^) 


*  This  lemma  is  proved  later. 
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which  satisfies  (2.20)  for  0  <  ^  <  ^''  <  ^'.   To  verify 
(2,21)  for  w   define  Q  =  (hJJ)  and  apply  Lemma  2*. 

co(r)lQn|r   [1  +   IFnglr-i  ]^^^^ 

>nlr  <  Zr 

(b+rm)^^ 


From  (2.18)   we  have 


J. 


l^nir 


1  <v 

<  2c(r)   +  Cin^  <     c(r) 


for     ^     small ,   0  <  ^  <  jl  <  n"  .     Thus 

Wnlr  ^   (2/b)'"'\5(r)co(r)   2^^^^    =K 

Repeating  this  procedure  we  generate  a  sequence  {w  }  which 
satisfies  (2.20)  and  (2.21)  for  Q  <  \i  <  '^. 

To  show  convergence  define  6u^^  (9)  =  u^^  (©ji-i) 
-  u^(e,^),   6G^^^  =  G^^^  -  G^  and  so  on.   Then  from  (2.22) 

6u^^^(e,)  -  (1-2mA)6^^^u(0)  =  ^[6G^^^  -  u^  6F^^J 

6v^^^(e.)  -S6v^^^(e)        =  6H^^^  -  v'^  6F^^^ 

Gi  =  e  +  T  +  ^Vi^®'^^^ 

Defining  Q   as  before,  we  have  from  the  first  estimate 
of  Lemma  2 

Using  (2.21),  (2.23)  and  (2.18)  we  have 


This  lemma  is  proved  later. 
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for  [1     sufficiently  small,  0<ti<|a  <ia,   i.e.,  uniform 

convergence  of  the  w  .   For  the  9  derivative  w 

n  n 

0  <  X  ^  r-1  v;e  have  from  Lemma  3  (with  m  =  1,  Xi  =9, 

-t  =  r)   and  (2.21) 

(O    (\)  I'ix/r) 

l^n-^p  -^n   lo^  2K  c(x,r)iv;^^p  -  w^]^ 

Hence   {w  }   converges  to  a  solution  y     of  (2.19)  and 
w(e,ia)  g  C"^"  (e)  n  C((j).   The  Llpschitz  condition  follows 
from  passing  to  the  limit  in 

Iw^^^^e',^)  -  v,^^^(e,^)|  <  Iw^l^+je'  -e|  <  Kje'-ei. 

Uniqueness  follows  by  applying  the  convergence  argument 
to  the  sequence   {wi ,W2 jWi ,W2 > . . .  }  where  wi   and  wa 
are  two  supposed  solutions.   Finally  we  note  that 
Iwl    <  K.  Thus  the  theorem  is  proved  with  a  =  u. 

Lemma  1.   (Solution  of  Linear  Functional  Equation) 
Consider  the  equation 


uOi)  -  (l-2iaA)  u(e)  =  nG(9) 


(2.2?) 


vOi)  -  sv(0)      =  HCe) 

©1  =  9  +  T  +  ^FC9,la)  (2.26) 

where     t     is   constant,      0  <  \i  <  \i' ,      A>0     and 
|Sv|    <  aolvj      for   some     Oq,      0  <  ao   <  1.     Define 
b   =  min(2A,    1-ao)      and     ra  =   inf^     F^    .      Assume 

1.  b  +  rm  >  0     for   some   integer     r  >  0 

2.  F,    G,   H   €   C^(9)    n   C(n)      and  have  period     2tt   in   9. 


-  3?  - 

Then  for 

1    1 

0  <  la  <  n"  =  min[ ,  ,  n',  1] 

2A  lm| 

(2.2^-6)  has  a  solution  u(e,vi),  v(e,n)  e  c^Ce)  n  C(n) 

having  period  2v     in  0. 

Remark.   Since  F  is  periodic  in  9,   m  will  be 
<  0  and  hence  if  1  is  true  for  some  integer  r  =  ro 
then  it  is  true  Tor  any  smaller  integer  r  <  ro . 

Proof.  Each  equation  has  the  form 


Y(ei)  -  D(|a)  Y(8>  =  PC6,h) 


with 


1D(^)Y1  <  (1  -  iib)|Y| 

where  D(n)   represents  either  S((a)   or  1  -  2|aA. 
Let  I  be  any  closed  subinterval  of   C0,^'').   We  may 
solve  (2.26)  for  the  negative  iterates  9   ,  9_2>  ..♦ 


(2.27) 
(2.28) 


where  9 


^,  =  9   +  T  +  iaF(9  ^,n)   since 
-n+i    -n  -n"^ 


d9i 
d9 


>  1  -I-  lam  >  0  for  la  €  I. 


We  allow  D  to  depend  on  9  also  and  write 

Y(9)  -  D(9_^,^)  Y(9_^)  =  P(9_^,n) 
D(9_^n)Y(9_^)  -  D(e_^,n)D(9_2,ia)YC9_2)  =  D(e_^  ,^a)P(©_2»^^^ 


Adding  and  passing  to  the  limit 

oo        j-1 

Y(0,u)    =J2     TT    ^^Q-k'^'^^^^-J'^^ 
d=l       k=l 
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where  n,  _^  =  1.   From  (2,28)  we  see  that  the  series 
converges  uniformly  for  all  9  and  for  la  6  I.   If  r  =  0 
we  are  finished.   If  r  =  1  we  differentiate  term  by  term 
using 


de 


and 


— 9-- 
d9  • 

d 
d9 


JZ2- 


d9 


-1 


d0 


-j+i 


d0 


•9. 


d 


^1  ^  mFq^ 

The  differentiated    series  converges  uniformly  for  the   same 
values   of     9     and     n     since     for      jy|    =  1, 

(2.29) 


|D(n)yl           1   -  lab 
<     — <  c.  <  1 


1  +  uFq         1  +  um 

where  c  depends  only  On  I,  h  and  m.  Hence 
Y(9,n)  €  0^(0)  n  C(^).   If  r  >  1  assume  for  s  <  r 
that  Y  €  C^"^(0)  n  C(^).   Then  differentiate  (2.2?) 
s-1  times  with  respect  to  0 


(9^  - 


D(|a)Y^^"^\9) 


(1  +  uF^)^ 


(2.30) 


ViS^Q. 


X\) 


where  P  contains  derivatives  Y  '^',   0  5  X  <  s-2 

considered  as  known  functions.  Thus  P  6  C^(9)  n  C{^i) . 

(s-i ) 

V7e  now  consider  this  as  an  equation  for  Y    '  and 

repeat  the  above  procedure  with  D(ij)[l  +  laFQ]^"^  in 
place  of  D(^).  Using  the  fact  that  for  x  >  -1, 
(l+x)   >  1  +  sx  we  see  that  the  corresponding  condition 
(2.29)  is 
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Diw.)y\  1  -  nb 

< <  Ci   <  1 


[1  +  ^FQ]^       1  +  nsm 
i^ence  the   solution  of   (2.30)    is   in     C^  (9)    n  C(^i) ,      i.e., 

Y(e,n)  6  c''(e)  n  cCn).  Q.e.d. 

Lemma  2.    "A-priori"   Estimates.      Consider   equation 
(2.25-6)    of  Lemma  1  and  assume  ever:/thing  in  the  state- 
ment   of  that   lemma,     letting     Q  =   C^)      we  then  have  the 
following  estimates    for   the  solution     v;  =   (^)      of    (2.25-6). 


2 

lvj|o   <  -    Nlo 
b 

2   IQU 
'    '  (b+m)2 

and   in  general  for     1  <  s   <  r,      s     an  integer, 

,6(s) 


wL  < 


co(s)lQ|3Cl  +    IFqU.,]' 


2  -  (b+sm)^""' 

where,   as  before,      m  =  infg     F^.      cq (s)  is  a  constant 

which  depends   only  on     s.      6(1)    =  0     and  6(n+l)    =   6(n)   +  n. 

Proof.      V/e  will  use  the  fact  that  for     x  >  -1     and 
n  >  0     an  integer      (1+x)      >  1  +  nx.  (2.3O) 

From  the   first   equation   of    (2.25) 

juOJl    <    (1-iab)    lulo   +  n|G|o 

for  all  9i    and   in  particular   at   the  point  where      lu|o      is 
assumed.     Thus      |u]o<b      |G|o. 
Similarly  for     v 

]vlo   <   (l"-b)|vlo   +    jHlo 
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-1 
and  therefore   |vjo  <  b   |H|o.  Adding  we  obtain  the  first 

estimate.  For  the  second  estimate  we  differentiate  to 

obtain 

Noting  the  restrictions  imposed  on  ^  in  Lemma  1, 

i:i+nm)lu^^^lo  <  (l"Mb)ju^^^lo  +  m|G1i 

and  similarly  for  v.  Thus 

(1 )  .1 


w 


jo  <  2(b+m)   |Q|i. 


2|Q|i 

Using  b+m  <  b  <  1  we  have   Iwli  <  .   To  prove  the 

"  (b+m) 2 
general  statement  for  s  >  2  assume  it  true  for  s-1  and 

differentiate  the  first  equation  in  (2.25)   s  times 

U^^^OJd+nFg)^  -  (l-2nA)u^^^e)  =  ia[G^^^  +  g^] 

where  g   is  a  linear  combination  of  terms  of  the  form 
s 

^(P)-|-Y  F^'i^  (2.31) 

i=l 

where     0  <   p,q  <   s-1     and      1  <  x.    <   s. 

At   the  point      Gi      where      Iu^^^(ei)|    =    |u        |o      we 
have,   using    (2,30),      1  +  lasm  <   (l+iaFg)^  and   therefore 

lu^^^lo(l+^sm)    <    (l-^b)|u^^^Io   +  n[|G|g   +    jg^lo] 
or 

b   +  sm 
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Using  the  estimate  for  s-1  and  the  fact  that 
b+sra<...<h+m<b<l,  we  obtain 

In  estimating  g  ,   since  the  maximum  number  of 
terms  in  the  product  in  (2.3I)  is  s-1,  we  choose 
6(s)  =  6 (s-1)  +  s-1.   Repeating  this  for  the  second 
equation  in  (2,25)  and  adding  we  obtain  the  desired  result. 

Q.E.D. 

We  now  establish  an  inequality  which  is  useful  in 
proving  the  convergence  of  the  derivatives  g    , 
0  <  X  <  "^j   of  a  sequence  whenever  the  original  sequence 
{g  }  converges  and  the  sequence   [g  •^ }   is  uniformly 

bounded. 

I 

Lemma  ^.   Suppose  f(x)  =  (fij..o,f^)  €  C  (x)   for 

all  X  =  (xi5.o.,x  ).  Then  for  0  <  x  <  l 

,  l"(xA)    \/l 

P^f  lo  1  oix,^)    |f  jo       jfj^ 

for  positive  integers  \     and  ^,   where  D^   is  any 
partial  derivative  of  order  x.   c  depends  on  m  and  n 
but  not  f. 

Proof.   We  define  the  norm 

^r 


[f]^  =  sup  |D  f^l 


X,  1 


where  the  supremum  is  taken  over  all  r   order  deriva- 
tives. We  introduce  the  smoothing  operator  (J.  Moser  [19]) 
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T,^  f(x)    =  J   ...   J  K^j(x-x')f(x')dXi'...dXj^        (2.33) 


valid  for  any     N  >  0     and  all     x.     The  kernel  is   given  by 

Kjj(x)    =  n\(Kxi)    ...   K(NXjjj) 

where     K(t)    6  C"     for  all     t     and     K(t)    =  0     for      |t|    >  1, 

and  T     ^      .  /l       ^°^     ^  =  ° 

t     K(t)dt   =\ 

'^  I  0,      0  <  k  <  s 

with  s  an  integer  to  be  fixed  later.   K  is  constructed 
[19]  from  a  C"  function  j6.      We  assume  that  the  mode  of 
construction  as  well  as  /6     are  fixed.   Therefore  K 
depends  only  on  t  and  s.  For  (2.33)  the  following 
inequalities  can  be  verified  for  a,  s,  r  >  0: 

CV^r+a  ^  Ci^'^'S)N^[f]^  (2.3^) 

[f-Tjjf]p  <  C2(s)   N^^Cf]^^^         (2.3?) 

where  Ci  and  c^  depend  on  m  but  not  on  N  and  f. 
Letting  N  =  ([f]  /[fjo)'^^'^,  we  use  (2.3?)  v/ith  r  =  \, 
s  =  ^- X    and  (2.3"+)  with  r  =  0,  a'  =  X   to  obtain 

The  final  result  follows  using 

c3(n)pPrio  <  [r]p  <  |flp 

where  C3   is  a  constant  depending  only  on  n.      Q.E.D. 


CHAPTER  III 
'  BIFURCATION  -  INVARIANT  SURFACE  METHOD 

1.  Introduction 

In  this  chapter  \ie   treat  the  bifurcation  problem 
by  reducing  it  to  a  problem  of  perturbation  of  an  invari- 
ant surface  of  a  vector  field. 

We  consider  the  differential  equation 

dx 

—  =  F(x,ia)  (3.1) 

dt 

where  x  and  F  are  real  n-vectors,  n  >  3>  and  n  is  a 
real  parameter.   Suppose  that  for  all  sufficiently  small 
u,  (3.1)  has  a  periodic  solution  x  =  >i/(t,|a)  which  is 
asymptotically  stable  for  la  <  0  in  the  sense  that  n-1 
of  the  Floquet  exponents  have  negative  real  parts,  i.e., 
n-1  of  the  Floquet  multipliers  lie  inside  the  unit  circle 
in  the  complex  plane.  As  \i     increases  through  zero  we 
assume  that  a  conjugate  pair  of  these  exponents  a(^)   a^d 
a(p)  crosses  the  imaginary  axis  and  enters  the  right  half 
plane  with 

Re  a(0)  =  0,     Re  a'CO)  >  0        (3.2) 

while  the  remaining  n-3  stay  in  the  left  half  plane. 
V/e  assume  that  FCx,|a)  has  p  >  5"  continuous 
derivatives  with  respect  to  x  and  \s     in  a  neighborhood 
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of  ^     and  in  that  neighborhood   |F|   <  Mq  where  differ- 
entiation is  with  respect  to  x  and  ^  and  Mo  is  constant, 

2.  -t^eduction  to  Normal  Form 

We  will  now  transform  equation  (3.1)  into  a  form 
suitable  for  discussing  the  behavior  of  solutions  near 
the  trajectory  iji. 

Following  NemytskiL.and  Stepanov  [20]  we  introduce 
the  arc  length  s  along  ^     as  the  new  Independent  vari- 
able and  yi,...,y„_^   as  new  dependent  variables.   The 
y.   are  coordinates  in  the  hyperplane;   S   normal  to  iji 
at  the  point   s.  Without  loss  of  generality  we  assume 
that  the  trajectory  ^     has  arc  length  2n,   i.e.,  tui  =  1. 
In  these  new  coordinates  (3.1)  takes  the  real  form 

y  =  A(s,n)y  +  i(s,y,|i)  (3-3) 

where  F(s,0,^)  =  0,  F  (s,0,^)  s  0  .  The  functions  are 
periodic  in  s  with  period  2rr  and   »   means  d/ds. 
Note:   If  we  let  So  be  the  hyperplane  S  discussed  at 
the  beginning  of  Chapter  II,  then  the  mapping  T:  yo   •*  y 
is  just  y(2T7,yo5ia)  where  y(s,yo,M)   is  the  solution  of 
(3.3)  which  satisfies  y(0,yo,[i)  =  yo . 

It  is  known  [21]  that  if  a(0)  /^  n  (i/2),  n  =  odd 
integer,  there  exists  a  real  periodic  matrix  Q(s,[j) 
having  period  2v     such  that  the  change  of  variables 

y  =  Q(s,^)u 
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reduces  (3.3)  to 


u  =  C(ia)u  +  G(s,u,u) 


(3A) 


where  G  has  period  2tt     in  s,   G(s,0,ia)  =  0, 
G  .Cs,0,^)  =  0  and  CC^i)   is  in  the  real  form 


C(u)  = 


Re  aCu) 
Im  a(M) 


Re  a(n) 


S(u) 


The  roots  of  C(|a)  are  the  n-1  non-zero  Floquet  exponents 
associated  with   i|i .   By  assumption  the  roots   oCO)   of 
S(0)   satisfy  Re  o(0)  <  oo  <  0  and  Re  a^u)   satisfies 
(3.2).  We  also  assume  that  a  linear  transformation  has  been 
performed  such  that   (Y,S(0)Y)  <  ao(Y,Y)  where 


If  we  let 


(3.^)  becomes 


Y  =  col(u3  , . . .  5U;^_,  ) 


z  = 

z  = 

Y  = 


U].  +  iU2 


Ui  "  IU2 


a(n)z  +  Gi (s,z,i,Y,n) 
S(n)Y  +  G2(s,z,z,Y,^) 


(3.5) 


where  G^   (complex)  and  G2   (real)  have  period  2tt  in 
s,   contain  no  linear  terras  in  z,z,   and   Y  and 
G^(s,0,0,0,n)  =  0.   The  functions   a(ki) ,   S(n)   and  G^  have 
p  >  5  continuous  derivatives  in  a  neighborhood  of 
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(z,Y,^)  =  0  which  we  assume  to  be   l(j|  <  u  , 
|2|^+|Y|^<2.   We  assume  the  G.   to  be  expanded  into 
Taylor  series  in  z,  z,  and  Y.   up  to  terms  of  degree  h 
plus  remainder.   As  in  Chapter  II  we  will  transform  (3.5') 
into  a  form  in  which  terms  of  weight  <  3  occur  in  a  par- 
ticular way.   We  recall  that  for  the  term  c(s,^)   z  z  Y."^ 
the  weight  t  is  defined  to  be  t  =  k  +  -C  +  2m.  This 
problem  of  removing  nonlinear  terms  from  a  differential 
equation  by  coordinate  transformation  has  been  treated  by 
many  authors.   CSee  G.  D.  flirkhoff  [22]  and  the  references 
given  in  Chapter  II  prior  to  Theorem  1.) 
Theorem  h.      In  (3.?)  assume  that 

1.  Re  a(0)  =  0,   Re  a'  (0)  >  0   ' 

2.  the  eigenvalues  of  S(0)   satisfy  Re  o(0)  <  0 

3.  a-(O)  =  V  i/3j  V  i/^5  V  an  integer 
^.  G.  6  C^  (s,z,z,Y,^) 

Then  there  exists  a  uo  <  u  depending  only  on  aCu)  and 
S(y)  such  that  for  \\j\  <  1^05  there  exists  a  transforma- 
tion 

z  =  cw  +  P(z,w,w,W,c,n) 

(3.6) 

Y  =  cV;  +  P(s,w,w,V/,c,|j) , 

c  constant,   0  <  c  <  1,  which  carries  (3.^)  into  the  form 

w  =  a()a)w  +  c2g(^)|w|^w  +  Ri^(s,w,w,W,c,|a) 

(3.7) 

W  =  SCn)W  +  R3(s,w,w,W,c,n) 


-  If?  - 

where  P  and  P  are  polynomials  in  w,v/  and  W  having 
no  terms  or  degree  zero  or  one,  and  the  R,   are  functions 
whose  Taylor  expansions  about   (w,w,W)  =  0   contain  no 
terms  having  weight   <  k.   The  mapping  is  one-to-one  in  a 
neighborhood   |wj^  +  |W|^  <  2  and  all  functions  of  s 
have  period  2n  in  s.  The  transformation  is  unique 
once  c  has  been  fixed.   In  case  3  is  violated,  the  fol- 
lowing is  true. 

If  aCO)  =  0  then  (3.7)  should  read 

^1   =  a(^)vJ  +  c[f  (^)w2  +  g(^)ww  +  h(ia)w2]  +  R3 

W  =  S(^i)W  +  R3  (3.8) 

where  R3  and  R3  are  functions  of   (s,w,w,W,c,^) .  ** 
If  a  Co)  =  iA  then  (3.7)  should  read 

w  =  a(n)v;  +  c2[p(^)(wj2w  +  Y(|a)e^^w3]  +  Ri, 

(3.9) 

w  =  s(^)v/  +  R3 

where  R3  and  R14.  are  functions  of   (s,w,w,V/,c,ij) . 

Proof.  Throughout,  subscripts  will  have  the  same 
significance  as  above,  i.e.,  R,   contains  no  terms  of 
weight  <  k.   We  write  the  first  equation  of  (3.5)  as 

2  =  a(^)z  +  U  (s,z,z,n)  +  .  .  .       (3.10) 


*  "Weight"  is  defined  in  Chapter  II,  Section  3. 


*tt 


For  this  assertion,  hypothesis  M-  may  be  replaced  by 
G^  €  C3. 
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where  U   contains  all  the  terms  of  degree  p, 

2  <  p  <  n,   in  z,z.   We  write  U   displaying  a  par- 

ticular  n   degree  term 

U^  =  U(s,z,z,n)  +  b(s,n)z  z'^,   k  +  ^  =  n 
where  oo 


bCs,^)  =   >   b   Cii)9 

m=-oo 

Define  dCs,n)  =  0  if   (k-1)  aCO)  +  I   aCO)  4   ni, 

n  an  integer,  and  d(s,^)  =  b^'^\n)e~^^^  if   (k-1)  a(0) 

**•  I   a(0)  =  ni.  Then  for  sufficiently  small  |j  the 

equation 

Y  +  [(k-l)a(^)  +  -ta(^)]Y  =  b(s,u)  -  d(s,^) 

has  a  periodic  solution  y^s,h)  having  period  2Tr  in  s. 
The  real  part  of  the  bracketed  term  is  /  0  for  ^  /  0 
and  therefore  YCs,ia)   is  uniquely  determined.  The  trans- 
formation 

z  =  w  +  y(s,^)w^w'^  (3.11) 

then  carries  (3. 10)  into 

w  =  a(u)w  +  U(s,w,w,n)  +  d(s,^)w  w  +  .  .  . 

where  the  dots   contain  no  terms   of  degree  <  n     in     w     and 
w.      If     a(0)  ^  V  i/3?    V   i/^     then  successive  applications 
of   (3.11)    for     n  =  2  and  3   leads  to  the  form 
w  =  a(u)v;  -^   g(u)|wl    w  +   .    .    .■    .      We  write  this   as 

w  =  a(M)w  +  e(n)|w|^w  +  Y«v(s,w,w,m)   "♦■  Gif  (3.12) 


'  h7  - 

■where     v  =  a(s,u)v}  +  b(s,n)\j,      a     and     b     vectors.      Con- 
sider the  vector  differential  equation 

Y  +    [S'(^)    •*■   {Ck-l)a(n)   +  -ta(^)}I]Y   =  hCsj^) 

where  prime  denotes  transpose,  I  is  the  identity  matrix 
and  h  is  periodic  in  s.  For  n  sufficiently  small 
this  equation  has  a  unique  periodic  solution  yCs,^)  and 
the  transformation 

w  =  C  •*■  Y  •  YCs,^)C^C'^  (3.13) 

carries  (3. 12)  Into 

C  =  aCu)C  +  pCn)|C|'C  +  Y  •  Cv(s,c,C,m)  -  hCs,u)c^c^]  +  Gi, 

applying  this  for  k  =  1,  t  =  0>  h  =  a(s,ia)  and  k  =  0, 
-t,  =  1,  h  =  t)(s,ij)  we  obtain  the  form 

After  carrying  out  the  above  transformations  on  the  second 
equation  of  (3.5')  we  write  that  equation  as 

Y  =  S(ia)Y  +  v(s,C,C»^)  +  H3         (3.1^) 

where  v  =  riCs,n)c^  +  ^'zQQ   +  ^'^V'     ^o^  M  sufficiently 
small  the  vector  differential  equation 

Y  +  CikaCn)  +  ^,a(^)3I  -  S(h)]y  =  h(s,|a) 

with  h  periodic  in  s,  has  a  unique  periodic  solution 
yCsj^j.)  and  the  transformation 


^-.l 
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Y  =  X  +  YCs,n)rC 
carries  (3.19)  into  the  form 

X  =  S(n)X  +  vCs,C,C,l^)  -  h(s,M)C^C^  +  H 


Applying  this  transformation  for  k  +  t  =  2  we  obtain 
the  form  X  =  S(n)X  +  H3 . 

The  previous  transformations  may  be  combined  into 

2  =  C  +  Q(s,C,C,X,n) 
Y  =  X  +  Q(s,c,C,X,n) 

which  is  one-to-one  in  a  small  neighborhood   |^|^  +  |X|^ 
<  2c,   0  <  c  <  1.   The  transformation  ^  =  cw,   X  =  cV/ 
then  gives  the  desired  result  (3.7). 

In  the  case  a(0)  =  0  we  apply  (3. 11)  only  for 
k  +  ^  =  2,  skip  (3.13)  and  proceed  from  (3.1^)  as  before. 
In  the  case  a(0)  =  i/^  the  procedure  goes  as  before 

except  that  in  equation  (3. 12)  there  now  appears  a  new 

/  ^  is-3 
term  yK\i)e      w  .  Q.E.D. 

3.  Non-resonance  Case:   Cs(Q)  /  v  i/^  1  v  iA-.   v  an 
inte.gjer 

In  this  section  we  obtain  the  results  of  Chapter  II 
by  the  method  of  perturbation  of  an  invariant  surface. 
The  problem  is  reduced  to  solving  a  certain  partial  dif- 
ferential equation  defined  on  a  torus.  The  actual  solu- 
tion of  this  equation  is  carried  out  in  Chapter  IV. 
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The  following  theorem  states  that  as  ^  increases 
through  zero  the  periodic  solution  ^  bifurcates  into  an 
asymptotically  stable  torus. 

Theorem  5.      In  (3.5')  assume 

1.  aCO)  /  V  i/3)  V  lAj  V  an  integer 

2.  Re  aCO)  =0,   A  =  Re  a'(0)  >  0  and  the  matrix 
S  =  S(0)   satisfies   (Y,Sy)  <  aoiY,Y)      for  ao  <  0  a 
constant 

3.  a(^),  S(n),  Gj^eC^CsjZjZjYjn),  I     an  integer  >  5 
h.   In  equation  (3.7)  B  =  Re  g(0)  <  0. 

Then  given  any  integer  r,  2  <  r  <  ^,  there  exists 

M-,  <  l-t   such  that  for  0  <  |a  <  n  ,   (3.?)  has  an  asymp- 
^  it  r 

totically  stable  invariant  torus 

z(s,e,ia)  =  aoj^   e   +  ur(s,e,n) 
TaCn):  (3.20) 

Y(s,e,(a)  =  ^g(s,e,^) 

where  ao  =  J-a/b,      f  and  g  are  defined  for  0  <  ^  <  ^^, 
have  period  2n  in  0  and  s, 

r,g€C^"^(e,s)  n  Lip^"^(9,s)  n  c(n) 

and   l^|p>i)   l^lr-i  -  ^^'^^      with  c(r)   a  constant  depend- 
ing only  on  r.   In  general  n  -  0  as  r  increases,  p 
depends  on  I-Iq   (defined  after  (3.2))  and  the  terms  of 
degree  <  3  in  the  expansion  of  (3.5)  about  (z,z,Y)  =  0. 
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Proof.   In  (3.6)  assume  c  to  be  fixed  such  that 
the  transformation  is  one-to-one  for   |v/|^  +  |W|^  <  2. 
Note  that  c  depends  only  on  terms  of  degree  <  3  in  the 
expansion  of  C3.?)  about   Cz,z,Y)  =0.  We  no  longer 
indicate  the  dependence  on  c  in  the  functions. 

We  restrict  attention  to  a  neighborhood 

w  =  az,     W  =  a^Y  (3.21) 

where   ]z|^  +  |Yj^  <  2  and  a,   0<a<l,isa  small 
parameter  to  be  fixed  later.   In  the  new  coordinates  (3.7) 
becomes 

z  =  a(|i)z  +  a^c^p(^)  |z|^z  +  fi  (s,z,z,Y,a,ia) 

Y  =  S(n)Y  +  f2(s,z,i,Y,a,n) 

with   |fi  li^  <  c(k)a3,  jfalk^^^^)^  where  differentiation 
is  with  respect  to  all  variables  except  a  and  la.  Note 
that  the  subscripts  on  the  functions  are  used  simply  for 
identification  and  have  nothing  to  do  with  "weight"  as  they 
did  in  the  previous  theorem. 

Expanding  a(^) »  S(^)   and  ^(n)  we  obtain 

z  =  Ca(0)  +  na'(O)  +  a2c2p(0)lz|2]z 
+  f 3(s,z,z,Y,a,n) 

Y  =  S(0)Y  +  fl,(s,z,z,Y,a,^) 

where 

Ifjl^  <  cCk)(a3+naH|a2)^  |fi.ll^<c(k)  (a+^) 
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Let     S(0)    =  S,    a(0)    =   iuj,   a'CO)    =  A  +  iA,    3(0)    =  B   +  IB 
and  choose 

0  <  l-i  <  lio    =  

-A 

Then  the  above  equation  becomes 

z  =  i[w  +  m/]z  +  nA[l  -  lz|^]z  +  r5(s,z,z,Y,n) 

Y  =  SY  +  f6(s,z,i,Y,p)  ^ 

BA 

where  id  -   A Izl^ 

B 
and 

1  i 

lf5lk  <  cCk)^iS       lf6|^  <  C(k)n 

The  unperturbed  equation  (fjjfs  s  0)  has  an  Invariant 
torus   |z|  =  1,  Y  =  0.   In  order  to  treat  the  perturbation 
terms  vje  introduce  coordinates  in  a  neighborhood  of  this 
surface  by  making  the  transformation 

z  =  (l+p)e^®,     |p|2  +  jYj^  <  -       (3.22) 

to  obtain 

p  =  -2nAp  +  [-iJLAp2  ^  r7(s,9,p,Y,u)] 

©  =  (ju  +  [^i(ci+C2p+C3p2)  +  f8(s,e,p,Y,u)] 

Y  =  SY  +  f9(s,9,p,Y,n) 

where  the  c.   are  constants.  V/e  now  treat  fg  and  the 
expressions  in  brackets  as  perturbation  and  note  that  the 
unperturbed  system  has  an  invariant  torus  Ta?   P  =  Oj 
Y  =  0,   which  is  asymptotically  stable  for  n  >  0.  The 
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problem  is  to  find  an  invariant  torus  t2  Tor  the  per- 
turbed system.  By  making  the  transformation 

p  =  Vu  u,     Y  =  v^  V  (3.23) 

we  obtain 

u  =  -2nAu  +  |a[G2(s,e)  +  G^  (s,0,u,v,M' )] 

&  =  oj  +  |i[F°  -I-  F(s,e,u,v,ia)]  (3.2^) 

V  =  Sv  +  G^(s,e)  +  G2(s,e,u,v,^) 

The  functions  on  the  right  have  period  2tt  in  9  and  s 

and  are  defined  for   1^1^+  l^j^  <  K^,   0  <  |a  <  fio 

=  min(^0)  IAk^)  where  K  is  an  arbitrary  constant  to 

be  fixed  later.   These  functions  satisfy  F°  =  constant 

and 

1G°I^  <  c(k);   |G^l^,    1F|^  <  cCk)u   . 

The  problem  of  finding  an  invariant  surface  for 
(3.2^)  is  equivalent  to  finding  a  solution  u  =  u(s,e,^), 
V  =  vCs,9,n)   of  the  partial  differential  equation 

su        au 

—  [uu  +  i^F]   +  —  +  2|aAu   =  uGi 
d9  ds 

3V  dv 

—  [ou   +  ijF]   +  —    -  Sv   =  G2 

where     Gj_  =  G?   •*"  G.      and     F  =  F°   +  F.      If  we  define 

^   =   C^'      ^=(^0'      D(w,9,s,^)    =      (^/^^     j){U-^uF)—  +— } 

d"      as 

and  P  =  pQ  _2J,  the  equation  takes  the  form 
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D(w,s,e,|a)w  +  Pw  =  G(w,s,e,n)         (3-25) 

which  is  a  special  case  of  equation  ih,h)    of  Chapter  IV 
with  ra  =  2,  xi  =  9,  X2  =  s,   A?  =  F°   and   A^  =  1.   It 
is  shown  in  that  chapter  that  for  a  given  integer  r, 
2  <  r  <  -t,  there  exist  values  K  and  ja  such  that  for 
0  <  ^  <  ^,  (3.2^)  has  a  unique  solution  Jzi   v  =   w(©,s,n) 
having  period  2rr  in  9  and  s, 

w€C^"'^(9,s)  n  Lip^"\e,s)  n  c(n) 

and   lvi|  _  <  K.  The  form  (3.20)  then  follows  by  applying 
the  transformations  (3.23),  (3.22),  (3.21)  and  finally 
the  transformation  (3.6)  of  the  previous  theorem. 

We  now  show  that   T2   is  asymptotically  stable. 
Let  u*(s,^),  v*(s,ia),  e*(s,n)  be  a  solution  of  (3.2if) 
which  is  sufficiently  near  T2   for  some  value  of  s. 
Write 

u*(s,^)  =  u(s,9*(s,m),u)  +  6u 
v*(s,m)  =  v(s,e**(s,n),n)  +  6v 

which  defines  the  variations  6u  and  6v,  where  u  and  v 
represent  ?2.   Putting  this  into  (3.2^-)  we  obtain 

6Ci  +  2nA6u  =  n(Gi^-UQF^)6u  +  laCGiy-UgF^)  6v 

(3.26) 
6t  -  S6V  =  (G2u-)aVgF^)6u  +  (G2v-MVqF^)  6V 

where  the  terms  in  parentheses  are  functions  obtained  by 
replacing  9  by  9*(s,ij)   and  the  derivatives)  F   etc., 
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arise  from  use  of  the  mean  value  theorem,  e.g., 

^u  ^  •Pu(s,e*,u(e*,s,^)  +  e6u,  v(e*,s,n)  +  G6v,)a) 

where  0  <  e  <  1.  From  the  form  of  F  and  G^  given 
just  hefore  *w3.25^)  we  see  that  their  derivatives  with 
respect  to  u  and  v  satisfy 

|F  |o,  ...  <  Con  J    Co  =  constant, 

uniformly  in  a  small  region  R:  )6u|^  +  |6vl^  <  M  contain- 
ing f2.  I'hus  (3.26)  has  the  form 

6u  =  [-2iaA  +  Pi]6u  +  P26V 
6v  =  Qi  6u  +  [S  +  Q2]6v 

where  P.   and  Q.   are  continuous  functions  of 
(6u,6v,s,^)  which  satisfy   |P^lo  <  t)^  ,  |Q^Io  <  b^i  , 
b  constant,  uniformly  in  s  for   6u,6v  in  R.   If  we 
define  U  =  \bu\^   +  n | 6v | ^   and  m  =  min(2A,  -ao)   then 
for  0  <  u  <  min[l,  (m/8b)^]   it  can  be  shown  that 

• 

U  <  -^(m/2)U  and  therefore  U  -  0  as   s  -  co .    Q.E.D. 

h.   Resonance  Case   a(0)  =  lA 

V/e  now  discuss  a  case  in  which  a  particular  low 
order  resonance  occurs.   Since  the  frequency  of  ^     has 
been  normalized   ((jui=l),   we  are  in  the  case  in  which 
0)2  =  lA-j  i.e.,  for  solutiors  of  the  linearized  (Floquet) 
equation  the  frequency  of  motion  along  the  periodic  solu- 
tion i|)  is  four  times  as  great  as  the  frequency  about   \j;. 
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In  this  resonance  case  we  see  that  the  normal  form  of  the 
mapping  (Theorem  h)    is  given  by  (3.9)  in  which  a  new 
parameter  Y^l-i)  now  appears.  The  relative  size  of  the 
parameters  a''(0),  p(0)   and  y{0)     determines  several 
cases.  We  discuss  one  possible  case  and  show  that  if 
bifurcates  into  a  pair  of  subharmonic  solutions  one  of 
■which  is  asymptotically  stable.  More  precisely 
Theorem  6.   In  (3.5")  assume 

1.  a(0)  =  iA 

2.  Re  a'(0)  >  0,   Im  a'(0)  =  0  and  S  =  S(0) 
satisfies   (Y,SY)  <  ao(Y)Y)   ^or  some  constant  oo  <  0 

3.  a(^),  S(^),  G^6C'^(s,z,z,Y,|a)   and  periodic  in 
s  with  period  2tt)   and  in  (3.^1)  below  assume 

i+.  Re  g(0)  <  0,  Im  b(0)  =  0,   and   0  <  |y(0)  |  <|p(0) 
Then  there  exists  a  ijo   depending  on  Mo   (defined 
after  (3,2))  and  the  terms  of  degree  <  3  in  the  expansion 
of  (3.?)  about   (z,z,Y)  •-  0  such  that  for  0  <  ^  <  lao  ? 
(3.?)  has  two  subharmonic  solutions.  Each  has  the  form 

z  =  Vm  20  e    +  laf  (s,^) 

(3.^0) 

Y  =  ^g(s,(a) 

where  zo   is  a  constant,   f  and  g  have  period  8tt  in 

1. 
s,   f,g  €  C^(s,n^),  i.e.,   C^  as  functions  of  s  and 

^'^,   and  \^\h->    \Z\^   <   constant  (independent  of  \s) .      One 

solution  is  asymptotically  stable  and  the  other  is  unstable 
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/  -a 'CO)     \^ 

For  the  stable  solution  zo  =  J   and  for  the 

\3(0)  +  1yC0)|/ 


unstable  solution  Zo  =  ' e  »*■ 


i      TT 

g(o)  -  |y(o)i, 

Proof .   We  apply  Theorem  h   to  obtain  the  canonical 
form 
w  =  a(M)w  +  c2[p(^)|vll^w  +  Y(|-i)e'^^v;^]  +  Ri^(s,w,vI,W,m) 

v;  =  s(^)w  +  R3Cs,w,w,w,n)  (3.^1) 

valid  for   |w|^  +  |W|^  <  2,   where  R,   contains  no  terms 
of  weight  <  k.   Letting  a'(0)  =  A,  c^g(O)  =  B  and 
c^yCO)  =  C  we  have  A  >  0,  B  <  0  and  by  replacing  s 
by  s  +  constant  we  may  assume  C  >  0,  Hence  0  <  C  <  -B. 
V/e  then  make  the  transformation 

w  =  J[sQy  W  =  nv  (3.^2) 

~hA  /B  +  C        \ 

where      \Q\^   +    Ivl^  <  — -^ —     and     0  <  [j  <  min ,    1    , 

B   +   C  \-2A  / 

to   obtain 

i  _  _ 

C   =  ~C   •*■   |ah(s,C,C)    -'-  A(s,C>C,v,n) 

(3.^3) 
V  =  Sv  +  r2(s,c,c,v,n) 

1  jV^  p  p  Q 

h(s,c,C)   =  AC  +  Bld^C  "*■  Ce^^^^ 


and 

,"2 


3  J. 


|ri  I  It   <  K,a2  ^       I  ^2  K   <  K^2  J      K  =  constant   >  0 
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In  this  system  we  consider  the  terms  ti     ^'^^     tz 
as  small  perturbations  and  discuss  first  the  system 
obtained  by  dropping  them.   V/e  will  show  that  the  unper- 
turbed system  has  a  pair  of  subharmonic  solutions  having 
the  stated  properties  and  that  this  remains  true  if  the 
perturbation  is  taken  sufficiently  small.   We  thus  intro- 
duce  a  new  parameter   e  by  writing  fi  -  ne  f i , 
f  2  =  e'^fa  where  \^ A^  <   K.  Then  for  e  =  i-i  we  have 
C3.^3)  and  e  =  0  gives  us  the  unperturbed  equation. 

Let   €  =  0  and  make  the  transformation 

s 

C  =  Tie^"^  (3.^?) 

to  obtain 

V  =  Sv 

16 
or  in  polar  coordinates,  t)  =  re 

r  =  ^[A  +  CB  +  C  cos  ^©)r2]r 

e  =  -^Cr^  sin  ^-9  (3-^7) 

V  =  Sv 

The  stationary  points  of  this  equation  are  (see  Fig.  8) 

where  n 

//  -A   \  in^ 


Tl.  =  V  e 


'n 


^n  =' 


B  +  C 


.B  -  C. 


rr    n 


e  ^'^    2 


v^  =  0,     n  =  0,1,2,3. 
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Fig.   8 


Ti  =  re    plane 


-A 


<ti:   r^  = 


B  +  C  cos  i+e 


A  =  if 

B  =  -  ?/2 

C  =  3/2 


r  consists  of  eight  such  arcs  joining  the  outer 
singularities.   Only  one  arc  is  shown. 
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These  stationary  points  correspond  to  a  pair  of  periodic 
solutions  of  equation  (3.^3)  with   e  =  0: 


The  linear  variational  equation  obtained  from  (3.^7)  has 

the  matrix 

'  Ji   0 
J  =■ 

0    S 


where 


Ji  =  I         ^A   I     ^^  ^n 


and 


Ji  =  1  CA  )     at  Q^  . 


Thus  poCs)   is  asymptotically  stable  and  qoCs)   is 
unstable.   It  is  easily  seen  that  po  and  qo  are  the 
only  periodic  solutions  (other  than  the  zero  solution)  of 
(3. 1+3)  for   e  =  0. 

We  now  treat  the  perturbed  equation.   For  ^a  =  e  =  0, 
C3.^3)  '^as  a  family  of  periodic  solutions  Q   =   QqB        ,v  =  0 
having  period  Stt  where  i^o   is  an  arbitrary  constant. 
According  to  standard  theory,  e.g.,  [H],  for  sufficiently 
small  u  and   e  there  exist  initial  conditions   CoC[a,e), 
Vo(u,e)   with   Co (0,0)  =  Co   and  Vo(0,0)  =  vo   and  a 
periodic  solution  of  (3.M-3)  satisfying  these  initial 
conditions  for  s  =  0  if  for  ^  =  e  =  0, 
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Ce^"^   -  I)vo   =  0 

877         t  t  t 

H(Co,Co)    H    i"       e"^*^  h(t,Coe^^,    Coe"^^)dt   =  0 

and 

.       dCH,  H) 

J  =  r—  ]^  0. 

s(Co,Co) 

Evaluating  the   integral,    we   obtain 

H(Co,Co)   =  8rrCACo+BjCoi^Co+CCo^) 

Hence  the   first  two  of  these  conditions  are   satisfied 
if      (,^°)      is   one  of  the  points     P„     or     Q_.      Evaluating  J, 


vje  obtain 


J   =    1HJ2    »    |H.|2 


=     [(C-B)2    -    (B+3C)2] 

\B+C/ 


at      P^      and  ,  2 

n  /8nA\ 

J  = [(C+B)2    -   (B-3C)2] 

Vb-c/ 

at     Q    .     Thus     J  /^  0     since     0  <  C  <    |B|      and   there   exist 
periodic   solutions      p(s,^,e)      and     q(s,u,e)      such  that 

p(s,0,0)    =  pCs,^,0)    =  poCs) 
q(s,0,0)    =  q(s,ia,0)    =  qoCs) 

These  solutions  exist  for  all  sufficiently  small  ^  and  e 
and  in  particular  for  ^  =  e.   If  we  assume  the  functions 
in  (3.^1)  are  smooth  as  functions  of   (s,w,w,W,ia)   then 
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•? 
(3.^3)  is  smooth  as  a  function  of   (s,c,C^v,^  ).  Hence 

p(s,^,|J)  =  Po(s)  +  u^p(s,^),  \^\\^.   <   constant,  where  the  ccristait 

is  independent  of  p,  and  similarly  for  q(s,n,u).  Using 

the  estimates  (3."+^),  it  follows  from  standard  arguments 

that  p(s,n,^)  and  q(s,p,^)  have  the  same  stability 

properties  as  poCs)   and  qo(s)   respectively.  The  form 

(3.^0)  follows  from  applying  (3.^2)  and  the  transformation 

of  Theorem  h. 

5,    Another  Resonance  Case 

In  previous  sections  we  considered  cases  in  which 
a  conjugate  pair  of  Floquet  exponents  crossed  the  imagi- 
nary axis  as  ^  increased  through  zero.   We  now  consider 
the  case  in  which  a  real  exponent  a(n)  enters  the  right 
half  plane  and  satisfies 

a(0)  =  0,     a'(0)  >  0  . 

Since  uua  =  0   (notation  of  Chapter  I)  there  is  no  tend- 
ency for  solutions  of  the  linearized  equation  to  oscillate 
about  the  periodic  solution  if.     Hence  in  this  case  we 
expect  that  ^     bifurcates  into  another  periodic  solution 
having  the  same  period. 

We  repeat  equation  (3.^) 

u  =  C(p)u  +  G(s,u,u)  (3.?0) 

where  G(s,0,^)  =  0  and  G^(s,0,n)  s  0  and  assume  that 
CCu)   is  in  the  real  form 
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C(n)  = 


a(|-i) 
0 


S(n)  J 


where  S(|a)   satisfies  the  same  conditions  as  before. 
Equation  (3.^0)  is  assumed  to  be  valid  for   lii|  <  n 
and   I  ul  ^  <  2  where  a   is  a  constant  >  0.   We  introduce 
the  notation  ?  =  u^ ,  U  =  col(u2,...,u  )  and  write  CS.^O) 


as 


I  =  a(ti)§  +  GiCs,5,U,n) 


U  =  S(n)U  +  G2(s,§,U,^) 


V/e  write  the  scalar  function  Gi  as 


(3.?1) 


GiCs,5,U,n)  =  3(s,n)§2  +  .  . 


and  define 


2n 


B  =  — 
2n 


3(s,0)ds 


Theorem  7.  Assume  in  (3.51)  that 

1.  a(0)  =  0,  a'(0)  >  0 

2.  CU,SCO)U)  <  aoCUjU)   for  some  constant   ao  <  0 

3.  B  /^  0 

^.  Gi,  G2€CMs,?,U,|i) 
Then  there  exists  a  yi,  ti  <  li  )   such  that  for 
0  <  n  <  (i,  (3.JI)  has  an  asymptotically  stable  periodic 
solution 

5  =  nco  +  n^fCsjii) 
U  =  la^gCs,^) 


(3.52) 
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where  Co  =  -(a'(0))/B;  r,g  e  C3(s,^);  lfl3,  \s\^   <   con- 
stant uniformly  in  n  and  f  and  g  have  period  2n  in 
s.  The  constant  \1     depends  on  the  Mo  defined  after 
(3.2)  and  the  terms  of  degree  <  2  in  the  expansion  of 
C3.?l)  about  U  =  0,  §  =  0. 

Proof.   In  C3.5'l)make  the  substitution 

5  =  ^x,     U  =  ijY 

/2a'(0)y  /  IB 

jx|2  +  IY|2  <  I 1  0  <  n  <  J2[ 

'  '     '  '     \  IB]  /  \2a'(0: 

and  expand  a  and  S  to  obtain 

X  =  |aFi(s,x,Y,n)  C3.?3) 

Y  =  SY  -4-  ^F2(s,x,Y,n) 

where  S  =  S(0)  and 

F^(s,x,Y,^)  =  a'(0)x  +  ^(.s,0)x^   +  ^F(s,x,Y,ia) 

where   |F]3  <  constant  uniformly  in  n.   For  n  =  0  (3.53) 
has  a  periodic  solution  po :   x  =  Xo ,  Y  =  0  where  xo   is 
an  arbitrary  constant.   It  is  known  that  for  sufficiently 
small  ;i,   (3»5'3)  will  have  a  periodic  solution  p  (s) 
of  period  2n  which  for  la  =  0  reduces  to  po ,   if  Xo 
satisfies  2 

H(xo)  =  r  Fi(s,xo,0,0)ds  =  0 

and  H'(xo)  /^   0.  Furthermore  the  solution  will  be  asymp- 
totically stable  if  H'Cxo)  <  0.   We  see  that 
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H(xo)  =  2n[a'(0)xo  +  Bxq^]   and  therefore  Xo  =  -  Ca'(0))/B 
and  H'(xo)  =  -  a'(0).  From  the  assumed  smoothness  of  the 
functions  in  (3.?3)  vje  have  p  (s)  =  Po  +  npCs,ji)  where 
jpl3  <  constant  uniformly  in  [i,  Q.E.D. 

6.  Resonance  Case   '  a(0)  =  i/^      (An  Example) 

We  now  discuss  an  example  of  this  resonance  case 
for  the  sake  of  completeness  and  to  point  out  the  rather 
violent  type  of  instability  which  develops  as  the  exponents 
aC|a)   and  oTilT  cross  the  imaginary  axis.  V/e  consider 
the  case  in  which  the  second  equation  of  (3»?)  is  absent. 
By  transformations  of  the  type  (3. 11)  we  may  put  the  first 
equation  of  (3.5)  into  a  standard  form  similar  to  those 
obtained  in  Theorem  h, 

w  =  a((a)w  +  c6(|a)e  w  +  R(w,w,s,iji)  ^ 

where  R  is  a  function  v;hose  Taylor  expansion  about 

1j6 
w  =  w  =  0  contains  no  quadratic  terms  and  w  =  pe   .   V/e 

consider  the  case  in  which  a(u)  =  i/3  '♦"  i-i)   c6  (f-i)  h  1  and 

R  =  -  jwl^w, 

.    /i     \       is-2 

w  =  (  -  +  lijw  ->•  e  w  -  Iwl^w  (3. 60) 

i  -2 
Letting  w  =  ze  3  v;e  obtain 

z  =  ^iz  +  i'  -  \z\'   z  ^3.61) 

i0 
In  polar  coordinates  z  =  re    this  equation  is 

r  =  (|ji  +  r  cos  39  -  r2)r 
9  =  -r  sin  39 
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There  are   seven  stationary  points, 

r   =  0 

2tt 

P^:      r  =  1  +  0(^),  e  =  n— 

n  3 

n  2n 

Q^:      r   =  pi  +  0(^2)  e  =  _  +  rj_ 

""  3         3 

n  =  0,1,2;   and  0(   )      is  meant   as     la  -  0.**     For     n  >  0,    the 

origin  is  an  unstable  node,      P       are   stable  nodes  and     Q 

are   saddle  points    (see  Fig.    9).      The  rays     r  >  0,    0  =  n(TT/3)> 

n  =  0,1,..., 5^     are  invariant   manifolds  and  the  disc 

A:      r  <   5^/^     is   invariant   as     s  -  +oo .      Thus     A     consists   of 

six  invariant   subsets 

?  IT  n 

A,  :      0<r<-,  k-<e<   Ck+D- 

^  "  H  3  3 

k  =  0,...,5.   Applying  the  theroem  of  Poincar^  and  Bendixson 

to  the   A,   we  see  that  the  points  P   and   Q   must  be 
k  -  n        n 

joined  by  arcs  which  are  solutions  of  C3.6I).  Let  r  :r  =  y(6,h) 
denote  the  closed  curve  formed  by  the  sum  of  these  arcs  and 
the  points  P   and  Q  .   Then  y  ^^^s  period  2tt/3  in  0. 
In  the  w,s  space  the  original  periodic  solution  ^     is  the 
line  w  =  0  with  points  2nn  identified.  The  z  plane 
rotates  in  such  a  manner  that  when  s  increases  by  2tt,   r 
sweeps  out  a  2-dimensional  torus 

S      .  y 

T  :   w  =  yU  -  -,n)e^'°  (3.62) 

^  3 

The  points  P   sweep  out  a  stable  periodic  solution 
p  =  p(s,^a)   of  period  6n  obtained  by  setting  j6   =  s/3 


*  See  Tltchmarsh  [2V]  for  definition  of   0(  ). 
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Fl(r.      7 


\ 
\ 

\ 
\ 
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in  (3.62).   Similarly  the  Q   sweep  out  an  unstable 

periodic  solution  q  =  q(s,n)   obtained  by  setting 

1^  -   s/3  +  n/3.   As  n  >  0  tends  to  zero,   q(s,|a)  -  0 

and  we  may  say  that  q  bifurcates  from  ijr  as  n  increases 

through  zero.  However  as  ^  -«  0,   p  approaches  a  limit- 

i- 
ing  position  p(s,0)  =  e  3.  Hence  in  the  z  plane  r  has 

a  limiting  position  which  consists  of  three  spikes 

2n 
0  <  r  <  1,   9  =  n-T-,  n  =  0,1,2  and  we  see  that   t   does 

not  bifurcate  from  %     in  the  sense  of  the  definition  of 

Chapter  I,  i.e.,  it  does  not  collapse  onto  ^     as  la  >  0 

tends  to  zero.   We  also  see  that  for  arbitrarily  small 

la  >  0  there  are  solutions  starting  arbitrarily  close  to  % 

which  run  out  a  distance  1  +  0(n)  to  the  orbit  p(s,^). 

This  is  in  contrast  to  the  non-resonance  case  of  Section  3 

where  nearby  solutions  may  run  out  only  a  distance  of  the 

order  J^     to  the  torus  (3,20).  Even  in  the  resonance 

case  aCO)  =  i/V  we  see,  by  the  same  reasoning  above,  that 

the  unperturbed  equation  (3.^6)  has  an  invariant  curve  r 

(Fig.  8).   Thus  equation  (3.^3) >  with  fi   and  f2   omitted, 

has  an  invariant  torus  r   lying  at  a  distance  of  the 

order  J\x     and  therefore  solutions  starting  near  ^     can 

run  out  only  that  far.  The  problem  of  the  existence  of 

T   for  the  perturbed  equation  is  not  treated  in  this 

dissertation. 

2.   Rate  of  Grovjth  of  Bifurcating  Manifold 

We  have  discussed  several  cases  in  which  the  periodic 
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solution  i}r  becomes  unstable  and  a  new  stable  manifold 
bifurcates  as  n  increases  through  zero.  Since  the  stable 
oscillations  of  a  system  are  the  only  ones  observable  after 
a  long  time,  it  is  of  interest  to  know  how  these  new  stable 
states  change  as  the  parameter  ^  changes.   In  particular, 
how  far  away  is  the  stable  motion  from  the  original  peri- 
odic motion  i{t?  Such  information  is  given  by  the  quantity 
P  =  p((a)  which  is  some  appropriate  measure  of  the  distance 
from  lif      to  the  stable  bifurcating  manifold.*  Of  course 
for  ^  <  0,   ,1;  is  stable  and   p(n)  =  0.  Let  S  be  the 
normal  hyperplane  to  the  orbit   i|(  at  some  point  of   i|i. 
In  the  case  in  which  one  real  root  crosses  the  origin 
(Section  5")  the  bifurcating  manifold  is  a  periodic  solution 
which  intersects  S  at  some  point  P.   The  distance  from  P 
to  ij;  is,  in  the  first  approximation,   p(^)  =  (constant )|a, 
(Fig.  10).   In  the  non-resonance  case  of  Section  3  the 
bifurcating  manifold  is  the  torus  (3. 20)  whose  intersection 
with  S  is,  in  the  first,  approximation,  a  circle  of  radius 
p(^)  =  (constant )7ia,   (Fig.  H) .   Also  in  the  resonance 
case  a(0)  =  i/M-  the  stable  bifurcating  manifold  is  the 
subharmonic  solution  (3.^0)  which  intersects  S  at  a  dis- 
tance p(^)  =  (constant )yia  in  the  first  approximation. 

Thus  we  see  that  in  the  last  two  cases  the  bifurca- 
tion is  more  violent  in  the  sense  that  near  the  bifurcation 


*See  Minorsky  [1]  for  a  discussion  of  bifurcation  and  the 
function  p. 
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point,  p'(n)  becomes  infinite. 

For  the  example  of  the  last  section   (aCO)  =  i/3), 
the  subharmonic  solution  p(s,n)   intersects  S  at  a 
distance  of   1  in  the  first  approximation.   Hence   pCp) 
is  a  discontinuous  function   (Fig.  12) 


p(n)  = 


li  <  0 

^1  >  0 
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Fig--  10 
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CHAPTER  IV 
INVARIANT  SURFACES 

1.  Introduction 

In  Theorem  5  of  Chapter  III  we  saw  that  a  closed 
invariant  surface  of  a  vector  field  was  defined  as  the 
solution  of  a  certain  partial  differential  equation  of 
first  order  involving  periodic  functions.   In  this  chapter 
we  treat  the  problem  of  existence,  uniqueness  and  smooth- 
ness of  the  solution  of  a  slightly  more  general  equation. 
Thus  we  consider  the  real  periodic  quasilinear  system  of 
first  order 

D(w,X5m)w  +  P(x,u)w  =  G(w,x■,^)  (^.0) 


wnere     x        .._,...,  ^^^ 


=   (xi  ,. ..  ,Xj^),   w  =  col(wi ,. . .  ,w^)    -   (y)     with 
u  =  colCui , . , . ,u^  )      and     v  =  colCvi , . . . ,v^   ). 

/Pi(x,n)  0 

P(x,ia)    = 

\         0  PaCxjM)/ 


and 


G(w,x,|i)    =      (^O    =  QOCx)    +  GCw,x,n) 
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where     P.      is  an     n.    x  n.     matrix,      I,     is  the     n.    x  n. 
>]  J  J  J  J  J 

identity,      c   =  l/ia  >  1     and 

m 

Ed 
a    (w,x,n) 

v-1  V 

with     a        scalar  functions.      The  functions     a    ,      P     and     G 
are  defined  for 

(w,x,^)  e  wXrX{0  <  u  <  w] 

with   V/:   |wl  <  constant  and  R;   0  £  x.  <  R.     a 

period  parallelogram. 

From  the  coefficients  of  d  define  the  vector 

a(w,x,Li)  =  (ai,...,a  )  and  assume  that 

m 

(1)  a(w,x,^)   =  uj  +  ^A(w,x,^)      if     ni   >  0 

(2)  a(w,x,ya)    =  A(w,x,|a)  if     n^    =  0 
where     u;      is   a   constant   vector   and 

A(w,x,ia)    =  A°(x)    +  A(w,x,n)  (^.1) 

The  case  ni  >  0  is  the  "degenerate"  case  mentioned  in 
Chapter  I  in  which  the  flow  against  the  unperturbed  surface 
vanishes  for  n  =  0  (see  U-.2  below). 

Equation  C^-.O)  expresses  the  fact  that  w(x,n)  =  (^) 
should  represent  an  invariant  surface  for  the  differential 
system 

u  =  -|aPi(x,fj)u  +  ^Gl(w,x,^t) 

V  =   -P2Cx,fa)v  +  G2Cw,x,n)  (^.2) 

X  =  a(w,x,^i) 
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In  Chapter  I  we  considered  the  following  differ- 
ential system; 

X  =  f Cx,y,M) 

(^.3) 
y  =  g(x,y,M) 

where  x  and  y  are  vectors,   [i  is  a  small  parameter 
and  the  functions  are  periodic  in  x.  We  assume  that  for 
^  =  0  C^.3)  has  an  invariant  surface  y  =  j6oix).        The 
problem  is  to  find  an  invariant  surface  y  =  /$(x,ij)  for 
u  >  0  such  that  /^(x,0)  =  j6oix).      If  we  assume  )6o(,x)    =  0 
then  (H.3)  can  be  put  into  the  form  (H.2)  with  the  u 
equation  absent   (ni=0)   by  setting  y  =  py 
-g^Cx,0,^)  =  P2(x,ia)   and  a(v,x,(a)  =  f(x,y,M).  Then  we 
see  that   A°(x)  =  f(x,0,0). 

The  basic  idea  of  the  method  of  solution  of  C^-.O) 
can  be  illustrated  v;ith  a  single  equation  in  one  variable 

a(u,x)u  +  b(x)u  =  f(u,x)  ih.h) 

where  x  is  a  periodic  variable  and  b(x)  >  p  >  0,   0  a 
constant.   The  equation  is  first  replaced  by  a  similar 
equation  with  G°°  coefficients  and  right  hand  side.   Then 
it  is  linearized  by  substituting  an  approximate  C"  solution 
into  the  functions  a(u,x)   and  f(ujx)   to  obtain  the 
C"  equation 

a(x)u  -^'b(x)u  =  f(x) 

This  equation  is  made  elliptic  by  adding  a  higher  deriva- 
tive: 


-  7^  - 
-eu   +  a(x)u  +  b(x)u  =  f(x)         (^.5) 

XX  X 

with   e  >  0. 

The  existence  of  a  C"  solution  of  this  equation 
is  then  treated  from  the  standpoint  of  the  classical 
theory  (Lemma  5').  We  then  obtain  an  "a-priori"  estimate 
for  the  solution  by  multiplying  (^.5")  tiy  u  to  obtain 

-cCiCu^)^  -  u2]  +  ■^a(x)(u2)^  +  b(x)u2  =  f(x)u 

At  the  maximum  of  u^ ,   (u^)  ^  <  0  and   (u^)   =  0.  Hence 

we  obtain,  independently  of   e, 

1 
|u|o  <  -  \T\o 

0 

This  estimate  and  similar  ones  for  higher  derivatives  are 
then  used  to  prove  convergence  of  the  approximations  to  a 
solution  of  (M-.H-).  The  parameter  e  =  e(A)  goes  to  zero 
as  n,   the  number  of  iterations,  goes  to  oo  . 

2.   notation- 

Let   A(x,u)  =  (A.       .  )  be  a  k~dimensional  array 

J 1  • » •  J^^ 

of  functions  defined  for  x  g  R  and  [x     in  some  interval. 

Define 

\  '  I 

k 


(A,B)  =  >    A.     ,  B.     . 
L )    Jl  •  •  •  Jj^   Jl  •  ••  -J 


I A|o=  max  ] A(x,^) | 
x€H 

1A|  =  (A, A) 

i 
(A,B)2  =  J  (A,B)dx     and     ||A|I  =  CA,A)2 

H 
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For     r     a  positive   integer   let 

d  d  5 


a"  = 


ax       sx  Bx^ 

'C-1  -0  2  I  J. 


1  <  ^.    <  m,     denote  any     r       order  partial  derivative  and 
let  m 

[r]  ^^i       1    - 

The  r   derivative  of  A  with  respect  to  x  in  an  r+k 
dimensional  array 

J  1  •  o  •   Jv 

For  r  =  1  or  2  we  write  A   or  A  .   Then 


lA^^V  =1]  l^'^l'  =1^  Ca^A,a^A) 

Cr]         [r] 

Finally  define 

[Al   =    max    lA^^^jo 
0  <  k  <  r 

Note  that  the  value  of  this  norm  is  a  function  of  the 
parameter  ^.   Any  estimate  involving   |A|   is  understood 
to  be  true  for  each  value  of  |a  under  consideration. 

3.  Main  Theorem 

We  will  now  prove  the  existence,  uniqueness  and 
smoothness  of  a  solution  of  (^-.O).   From  the  vector 
A°(x)  =  (Ai,...,A°)  define  the  symmetric  matrix 


-  76  - 


aA. 


BX 


,A  0 


dx. 


Let   X^Cx)  denote  an  eigenvalue  of  this  matrix  and  define 

X   =   min  xAx) 
k,x  ^ 

Since  A° (x)   is  periodic  the  functions  sA?/ax.  have  zero 
mean  value,  i.e. , 

J  — ^  dx  =  0. 


R  dx. 


In  particular,  the  trace  of  the  above  matrix  has  zero  mean 
value  and  thus  for  some  k  and  x,  >^vCx)  <  0.  Therefore 
X  <  0. 

Theorem  3.   In  (^.0)  assume  that  for  all  x  e  R  and 

H5   Q  ^  ki  <  no, 

1.  (w,P(x,u)v)  >  9(w,w)   for  some  constant   p  >  0 

2.  3  +  rx  >  0  for  some  integer  r  >  2 

3.  Cvj,Pi  \})    "•  -^   div  A°  (x)  >  3o  >  0  for  some  constant 
Bo   and  all   |wl  =  1.   If  ni  =  0,   replace  Pi   by  P2. 

h.   PCx,^),  G(v;,x,^),  ACwjX,ia)  6  C'^(w,x)  fl  C(n)   for 
-t  an  integer  >  r  and   |w|  <  K,  K  sufficiently  large 


5.  maxi|A|^,  IGj^}  =  c(r,K,n) 

where  c  -  0  as  la  -  0. 

Define  M=max{jA°|^,  |G° j^,  sup^ jP |^} . 


'  IT  -■ 

Then  there  exist   constants  K  and  |a,   m  <  no   such 
that  for  0  <  n  <  la,   (^.0)  has  a  unique  solution 
w(x,m)  €  C  ~  (x,|a)  n  Lip^"  (x,n)  n  C(ia)   and   jw]^.^  <  K. 
K  depends  only  on  M,   g  +  r\  and  r.  \x     depends  on  r, 
^5  c,  3,  00   and  K  and  In  general  \i     decreases  as  r 
increases. 

Remark.   If  A°  (x)  -=  constant  as  in  Theorem  5"  of 
Chapter  III,  then  X  =  0  and  the  smoothness  of  the 
solution  of  (4.0)  is  limited  only  by  the  smoothness  of  the 
equation  and  the  smallness  of  \x.      Thus  if  ^  =  oo  in  (4) 
above  then  there  exist  intervals  0  <  (i  <  u  ,   n    <  ^ 
with  r  arbitrarily  large,  on  -which  the  assertion  holds  , 
i.e.,   v;  is  smoother  for  smaller  values  of  \x.      Reversing 
this  we  see  that  as  \x     increases,  the  number  of  deriva- 
tives which  the  theorem  guarantees  is  diminished. 

Proof.  Define   g '  =  3  +  rX   and  let  K  be  a 
constant 


2  \        ^ 


K>^-J   3M 


(3M,  2M)  +  1 


where  ^   is  the  polynomial  with  positive  coefficients 
defined  in  Lemma  4  (proved  later).   We  construct  successive 
approximations  w  (x,n)  with  v;o  =   0.   This  is  done 
inductively. 

Suppose  Wv,.!  (x>m)   satisfies 

w„  ,  €  ^^"(x)   for  each  fixed  n,  0  <  ^  <  u'  (4.30) 

w^.!?  o'^'^(x)  n  c(n)  (4.31) 


l"n-xlr<K  CH.32) 
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On  the  set   i|w|  <  K]  X  R  X  {0  <  |i  <  ^o }   approximate  the 
functions   A  =  A(w,x,n),  G  =  G(v;,x,n)   and  P  =  P(x,n) 
by  functions  A,  G,  P  e  C^CwjXjia)   and  consider  the  linear 
partial  differential  equation 

5^Cx,^)w  +  P(x,n)v;  =  G^(x,^) 


where 


G^(x,n)  =  G(w^_^(x,n),x,n) 


and  (^.33) 

1^^(X5^)  =  D(w^_^^'(x,|a),x,|a)  , 

i.e.,  in  the  expression  for  D,   the  function  A  is 
replaced  by 

A^(x,n)  =  A(v^__^  (x,ia),x,ia). 

We  assume  the  approximation  to  be   so  close  that 

lA      -A    I      +1G      ~GI+|P-Pi<M„        (^.3^) 
In         n'r        'n         nir        i  'r-n 

where 

1    /  3'    3o    3 
M  =  -  mini— —  ,  — —    ,  -— z  »  M 
^   n^    \12r   Vm   V''^ 

and  the  subscript  r  denotes  differentiation  with  respect 
to  X.   A   and  G   are  defined  by  (^.33)  with  the  tildes 
omitted.   V/e  modify  this  equation  still  further  to  obtain 
an  elliptic  equation  for  which  the  existence  and  smoothness 
of  a  solution  is  well  knov/n.   Thus  consider 


L^w  =  -e^ZiW  +  D^(x,^)w  +  P(x,|a)w  =  G^Cx,^)  (^.3?) 
where  a  is  the  Laplacian  in  m  variables,  e  -  0  as 
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n  -  00    and      e^  >   e^+^    >  0. 

Define     a     =  w  +  uA       and     use    (M-.3^)    to   obtain 
n  n 

Cw,    [P   -  idiv  a^C]w)    >  ^  min[3o ,3] («,«) 

for     n     sufficiently  small.      Hence  by  Lemma    5,    (^.3^)   has 
a    solution  w  =  ■w^(x,^)    €  C"(x)      for   each  fixed     ^.      If  we 
define     p  =  3   -  3 '/12     then  using    (M-.3^)   v/e  obtain 
(w,Pw)    >  p(w,vj).      Let      X   =  x(x,n)      be   obtained  from   X^(x,|a) 

just   as     X     was  obtained  from     A^  (x) .     Then  using    (^-.3^) 
and  the   fact   that      r|\-xl    <  r|A°-A^|-i_   <  rM^  +  2Krc(l,K,n) 
we  have     p  +  r\  >   3 '/2     for     \i     sufficiently  small.      V/e 
see  here   that   for   larger     r     we  need  smaller     \i.     Thus 
applying  Lemma  h, 


2\r+i 


l^nlr  ^ 


^nlrTT   MV^..    '    l^lr^ 


,3'/       '  ''  "■      II         I  '^xx 


r-2 


We  have 


l^nlr  ^  \  ■"    l^°lr   "    l^lr  <   2M  +  c '    <  3M 


A 
^xx 


r-2 


<  M     +    I  A°  I      +    I  A I      <  2M  +  c '    <  3M 
-     n        '       Ir        I    ir  -  


and 


IP  I      <  M     +    IP!      <  2M 
I    Ir  -     n        I     ir  - 


for     u     sufficiently  small.      Differentiation  is   with  respect 
to     X     only,      G     and      A     are   as  defined   in    (M-.l)   with 
w        Cx,^)      substituted   for      w.      c        depends    on     r,   K     and 
c     defined  in  hypothesis   5  and   is   obtained  by  applying   the 
chain  rule  to     G     and     A;c'-0     as     ^i-»0.     Hence  from 
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the  way  K  was  defined,   |w  j^  <  K.  We  now  verify  that 
w   satisfies  (M-.31).  For  a  function  f^Cx,ia)  define 

6f„  =  i'„o  -  ^ni  =  ^J^y   1^2)  -  r  Cx,m).  From  (^.35) 
n    n2    ni    n  n 

L   fiw  =  6G  -  C6D^)w„,  -  C6P)w^T 
n2  n     n      n  nl       n.l 

where  L    is  L   with  ^i  =  H2«  From  the  estimate  of 
Lemma  ^, 

Since   |w  |   <  K  this  can  be  made  small  for  (J2  ~  m  small. 

From  Lemma  3 5  Chapter  II,  ■ 

1-  - 
l^Wj^'^'^lo  <  2Kc(j,r)|6w^|o  ^ 

for  0  <  j  <  r-1.  Thus  w^   is  continuous  in  n 

uniformly  in  x.   Since  w  €CCx)   for  fixed  )j,  we  have 

w  6C^~^(x)  n  C(|a).   Thus  for   ^'   sufficiently  small  and 

fixed  we  have  a  sequence   {w^}   satisfying  (U-.30)  -  CM-.32). 

To  show  convergence  define  6_.,  =  w  , .  -  w^.  Then 

n+i  n+-^  n 

from  (M-.35)     ^n+i^n+i    ~  ^     w^e^e 

«  =  K-^r^n   -  ^^n^^n-S-M^    "    ^  Vx "  ^n^^n 
From  the  estimate  of  Lemma  h 

iVxio^;  i»io 

We  have 


<  2KM     +  KiA    |o    1 6    1 
-         n  '    wl"    I    n 
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where  A  =  A(w   jX,|a).  Similarly 
n      n-i 

IG  ^   -  G  1   <  2M  +  |G  I   16  I 
and 

Thus  for     la     sufficiently  small 

1   n+j.  lo  -  '^  1    niQ  n 

wher  e 

.^n  =  ^^P"'\  "■  2p"'^<^^n-^n-^x^- 


Since     z"  q     <  oo    v/e  see   that      E**    1 6  ^    I      <  oo    and  hence 
1    ^n  0    '   n+i ' 0 

the     w       converge  uniformly  in     x     and     (a. 

i'etting     2  =  w  +   /-  w   ,      p  >  0     an  integer,  we  have, 
using  Lemma  3 5  a 

|z^^^|q<  2Kc(j,r)lzlo      ^ 

for  0  <  j  <  r-1.  Thus  the  sequence  {w  }   converges  to  a 
limit  function  w  6  C   (x)  (1  G(n).   The  Lipschitz  condition 

follov;s  from 

lWn^''"'^(xx,u)  -  w^^^"'^X2,fi)|  <  |v;^l^  |xi  -  X2I 

<  K  jxi  -  X2I 

and  passing  to  the  limit.  For  uniqueness  assume  wi  and 
W2  €  C^Cx)  n  C(n)  are  solutions  of  ih.h)   and  satisfy 


WjU  S   K.  Then  6w  =  W2  ~  Wi   satisfies  the  linear  equa- 
tion 
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D26W  +  P6w  =  6G  -  C6D)wi 
where 

D2  =  DCw2,x,h),  6G  =  G(vi2>x,^)  -  GCwi,x,h) 
and 

6D  =  DCv;2»x,ii)  -  D(wl,x,^). 

Now  apply  Lemma  h   in  the  c  =  r  =  0  case  to  obtain 

2 

-  |6w|o  <  |6G|o  +  |6Alo  K 

3 

=  |6Glo  +  |6AIo  K 


<  (|G  I   +  lA  1  K) j6wl 


1 

<  -  j  6w|o 

for  n  sufficiently  small  (Hypothesis  5) •  Thus   |6w|o  =  0. 

Q.E.D. 

h.   "A-priori"  Estimates  and  Solution  of  Linear  Problem 

We  now  derive  some  estimates  appropriate  for  a 
solution  of  the  real  linear  elliptic  partial  differential 
equation 

Lw  5  -eAw  +  DCx,ia)w  +  PCx,|a)w  =  G(x,|a)     ih.hO) 

where  x  =  (xi , . . .  ,Xjjj) ,   w  and  G  are  column  n  -  vectors 

with  w  =  Cy),  u  =  colCui , . . . ,u^  ),  V  =  colCvi , . . . ,v^  ), 

1 
/  cli   0  \      c  =  -  >  1 

D  =  Cd  =  (  1*^'        ^ 

^0    lo/      I.  =  n^  X  n.  identity 

0    d    J 

with  m 

d  =  y^  a  (x,^Jl) .      a   scalar  and 

v-i  V 
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Pl(x,^) 


P(x,m)  = 


0 


P2(x,m) 


,  P  .   an  n .  X  n .   matrix. 


From  the  coefficients  in  d  vje  have  the  vector 
a(x,  )  =  (ai,...,a  )   where   (1)   a(x,ia)  =  u;  +  ^A(x,^)   if 
njL  >  0  or  (2)   a(x,^)  =  A(x,n)   if  ni  =  0  where  uu  is  a 
constant  vector.  We  assume  A  =  A(x,^),   P  =  P(x,n)  and 
G  =  G(x,u)   are  defined  for   (x,n)  6  R  X  I-  where  R  is  a 
period  parallelogram  and  I-:   0  <  n  <  |a. 
Consider  the  symmetric  matrix 


1. 
z 


SA.    SA  ■ 


dx, 


ax. 


and  let  Xv^^>l-i) 
Define 


0   -  ^i-i 
k  =  l,...,m  denote  the  eigenvalues. 


\ '   =   inf  min  x  (x,^) 
H  k,x  ^ 

and  note  that   X '  <  0  due  to  the  periodicity  of  A  and 
the  same  argument  as  given  before  Theorem  8. 

Lemma  h.        "A-priori''  Estimates 

Assume   e  >  0  and  for  r  >  0   an  integer,  assume 

1.  A,  P,  G  e  C^(x)   for  each  fixed  ^  €  I-  and 

w  ^  C^Cx);  I  =   r+2  if  e  >  0,  -6  =  r+1  of  e  =  0. 

2.  (w,Pw)  >  0(w,w)   for   Cxjij)  €  R  XI-  where   0  >  0 
is  a  constant 

3.  0  +  rx'  >  0 

Then  2 ] G | o 

lw|o  <  


-  Qh   - 


and  in  general  for     0  <  s  <  r 

s 

IGI   ~~ 


MAxxi..2'   ip:^ 


for  all  n  6  I-  ,  where  for  k  <  0  |  j,  is  taken  to  be 
zero,  -u(a,b,)  is  a  polynomial  in  a  and  b  of  degree  s 
having  positive  coefficients  which  depend  only  on  s.  The 
estimates  are  independent  of  e. 

Proof.   Due  to  the  form  of  P  and  the  fact  that   A 
and  D  are  diagonal  operators,  (^.HO)  has  the  form 

-eAu  +  cdu  +  PiU  =  Gi  (^.^2) 

-eAv  +  dv  +  P2V  =  G2  (^."+3) 

We  have   jw|o  <  |u|o  "*"  Iv|o  <  2|w|o5 

|P|o  <  |Pilo  +  |P2|o  <2|P|o   and   (v,P.v)  >  0(v,v)   for 
3  =  1  and   2,   We  obtain  the  estimate  for  (M-.^3)  and  note 
that  the  estimate  for  (,h,h2)   follows  by  replacing  a  by 
ca. 

Multiply  (^.M-3)  by  v   (scalar  product)  and  use  the 
identities 

(v,Av)  =  iA(v,v)  ~  jv^  f 
(v,dv)  =  -^dCvjV) 
to  obtain 

-e[iA(v,v)  -  [v^I^]  +id(v,v)  +   (v,P2v)  =  (v,G2). 

Let   xo   be  a  point  at  which   (v,v)   is  a  maximum,  i.e., 
at  Xo  ,  (v,v)  =  \v\§.      Since  a(v,v)  <  0  and  d(v,v)  =  0 


-  8^  -    ^ 

at  xo  we  obtain  (vjPgv) 1  _^  <  (v,G2)l      and  therefore 

elvjo  <  |G2|o. 

To  obtain  the  general  estimate  (^,^1)  we  assume  the 
solution  has  a  sufficient  number  of  derivatives  and  apply 
the  same  argument  to  the  differentiated  equation.  The  rule 
for  applying  d^  to  a  product  fg  of  two  functions  is 


d^(fg)  =  f  Cd%)  +  Y^ 


s 

3f    /  s-i 


where 


j=i  ^% 


3j  E)    * 


s-1 


ax.      ax,   ax        ax 


and  the  dots  (present  only  if  s  >  2)  represent  a  sum  which 
is  bilinear  in  f  and  g  involving  derivatives  of  f  of 
order  >  2,  Thus 

a^Cdv)  =  d(a^v)  +  Q^  H-  R^ 

where  s   m  / 

r— T  IT— -1  aa   a   /  s-^ 

k  k  '\^  '\  ^ ' 

and 


R  =  )      [v^P^^  a  ^^^] 

S    / t  V    -^ 

V 

p+q=s+l 
l<p<s~l 


where  R  =  0  for  s  <  1,   Q  =0  for  s  =  0  and  the  ^ 
s  --   '   >s 

notation  denotes  a  linear  combination  of  terms  a^va^a  .  Now 

V 

apply  a   to  (^.M-3),  multiply  by  a"'v  and  sum  to  obtain 
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where 


H  =  >   (a^Vjd^Ga  +  Rg) 
[s] 


with 


R 


-R  + 
s 


}    ;  [P.^^'  »  v«p'] 


t+p=s 
0<p<s-l 


We  now  show  that 


Y^   <;s^v,P2    S^v  +  Qg)    >   (p+sx')    |v^^^|2    ih.h5) 

[s] 


From  the  assumption    (2) , 


\  (^''v,   Paa^v)    >   3   )        (S^v,d^v) 


[s] 


[s] 
3    \^^^'\ 


CS),2 


For   the  second  part  ,  S  ,  we  have 


n2 


n2      s 


S   H  2_^   (ci"v,Q^)    = 

[s]  i-1    [s] 

ra  m 


:{> 


-♦    aa 


_iL 


bx 


p-l,.*.,s 


ax 


^%'«s^l 


^s-1 


V 


s 
S   v.l 


The   term  in  brackets   is   a  quadratic   form, 

m 
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vjiiere 


aa 


a 


_i^ 


Vjt 


j         dx 


I. 


_   "  s-i 

V 


dx 


.^r^'i 


^j 


and     e     is  the.  vector  of  subscripts 


6    -     vi,     ^1)     •••     't'-;-i''^-j+j_» 

From  the  way  in  which  \'      is  defined, 


•  •   •       ^  c  * 

5 


m 


K  >   \' 


ke 


and  therefore 


s        n2 


m 


S   >  \ 


=   \' 


j=l  i=l     k^^p=l  ^""k 

s  P  ^  -^ 


■which  proves    (H.>+5).      We  now  repeat   the   argument   after 
(4. "+3)   where     xo      is   a  point  at   which 


[sj 


is  a  maximum  and  obtain 

Using  the  fact  that 

[3  +  sx']  <  [3+  Cs-1)\']  <  ...  <  3  +  X'  <  p  <  lP|o   (^.^6) 
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v;e  write  this  as 

Cp  +  sx')^''"|v^^^j^<  p^  IG|^  +  [p  +  (s-l)x']  IR3I0 

We  now  prove  (M-.M-l)  inductively.   Assume  it  is  true  with   s 

replaced  by  0,  1,  ...,  s-1.  The  expression  R   is  linear 

in  V  and  each  term  contains  a  factor  d  v  with  p  <  s-1. 

Using  the  induction  assumption  and  C^-.M-S)  we  obtain 

s  ^ 

where  -jj-'      is  a  polynomial  of  the  type  described  after 
ih.hl).      (Note:   [a^]  =  \\sA^\    <    \A^\      since   0  <  ^  <  1. 
In  deriving  the  estimate  for  u  in  ("+.^2)  we  get 
lea  I  =  ]A  ]   and  hence  the  estimates  are  exactly  the  same.) 


From  (^-.M-l),  with  s  replaced  by  s-1,  we  have 

s 


vjhere  we  have  again  used  (^-0^6)  and  the  fact  that 
1  Ip  ^-  1  Iq   for   p  <  q. 

Now  define  -n     to  be  the  polynomial  whose  coeffici- 
ents are  twice  the  maximum  of  the  corresponding  coeffici- 
ents in  4'  and  u"  •      I'hen  O.^fl)  follows  since 

IVI3  =max(jv^^^|o,   JVI3.,). 

In  regard  to  the  assumption  made  on  i     in  (1), 
note  that  in  (M-.H-M-)  for  s'=  r  we  need  v  e   C  if 

e  >  0  and  v  €  C^"*"^   if   e  =  0.  Q.E.D. 
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Finally  comes,  the  problem  of  existence  of  a  solution 
of  equation  (^-.i+O) ,   Lw  =  G.   v;e  no  longer  employ  the 
maximum  norm  used  in  the  solution  of  the  quasi-linear 
equation,  but  introduce  the  Hilbert  space  of  square  inte- 
grable  functions  having  period  parallelogram  R.  By 
obtaining  an  "a-priori"  estimate  for  a  solution  of  the  adjoint 
equation  we  prove  the  existence  of  a  "weak"  solution  of 
Lw  =  G.  The  proof • follows  that  of  Friedrichs  [12;  p.  3?3]. 
Certain  difficulties,  which  arise  in  considering  domains 
with  boundaries,  disappear  due  to  the  periodicity.   The 
C"  differentiability  of  the  weak  solution  then  follows 
from  a  result  of  P.  D.  Lax  [23].  The  following  lemma  is 
a  statement  of  these  classical  results. 

Lemma  5.    ~      (Existence)  In  (^,^0)  assume 

1.  e  >  0 

2.  (w,  [P(x,^)  -  i  dlv  a(x,n)C]w)  >  p(w,w)   for 
(x,n)  6  R  X  I-  and  p  >  0  a  constant 

3.  for  each  n  e   I~,   A,  P,  G  €  C^Cx) 

Then  for  each  fixed   n  €  I_  ,  (^+.^-0)  has  a  unique 
solution  w(x,|j)  £  C'°(x).. 

Proof.  For  the  equation  Lw  =  G  we  have 

r  XwjLw)  =  r  (w,G) 
R         R 

and  integrating  by  parts 

sf  |w^l^  ■»■  f  (w,[P  -  i  div  a  G]w)  =  r  (w,G) 
R         R  R 
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where  periodicity  was  used  to  drop  out  the  boundary  terms. 
From  this  we  obtain 

■     llwll  <  -   IIGII  CM-.^7) 

P 

For  the  adjoint  operator  LV  ='  ~  eA^  ~  D(x,|a)/  +  P(x,n)/ 
with  P  =  P  -  Cdiv  a)C,  we  obtain  the  same  estimate 

W.  ^  -    ||L*;2il|  (^.1+8) 

P 

Let  LsCR)  denote  the  Hilbert  space  of  functions  f (x) , 
X  e  R,  such  that   l|f  |j  <  oo  .  ■  Using  C^.M-8) ,  for  each  fixed 
H,  we  prove  the  existence  of  a  "weak"  solution  of 
Lw  ■=  G,  i.e.,  a  function  w(x,^)  €  LgCR)   such  that 
(w,  LV)2  =  (G)/)2  ^or  all  /  in  a  dense  subset  of 
LzCR).  For  define, 

Q  =  [q:   q  =  L*/$,  j6   €  C~Cx)  ,    x  6  R,   /^ 

periodic  with  period  parallelogram  R} 
On  Q  define  the  functional  ii;Cq)  =  (G,/)2.   The  estimate 
(H.^-B)  implies  uniqueness  for  the  adjoint  equation  and 
hence  ^     is  v;q11  defined,  i.e.,  given  q  €  Q>   there  is 
exactly  one  /6   e  C"(x)   such  that  q  =  LV.   i(i  is  bounded 
and'  linear  due  to  C^.^8)  and  the  linearity  of  L^.   Thus 
by  Riesz's  representation  theorem  ijiCq)  =  (w,q)2  ^ox   some 
w(x,^)  £   LsCR),   i.e.,   (G,/$)2  =  (w,LV)2. 

In  the  previous  reference  to  P.  D.  Lax,  it  is  shown 
that  for  such  a  weak  solution  wCx,^),  if  G  and  the  co- 
efficients of  L  have  derivatives  in  x  up  to  order  t 
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which  are  in  L2(R)  then  v/(x,n)  has  derivatives  in  x 
up  to  order  t+2  vfhich  are  in  LgCR).  This  implies,  by 
Sobolev's  lemma,  that  w(x,|a)  €  G'^^Cx)  where  r  =  t  -  [nT/2] 
+  1.  Hence  in  our  case  w(x,|a)  e  C^Cx)   for  each  fixed  n. 
Uniqueness  follows  from  the  estimate  (^A7). 

Q.E.D. 

5,   An  Example  of  Non-parallelizable.  Flow 

The  following  is  an  example  of  the  situation 
described  in  Chapter  I  in  which  the  flow  on  the  unperturbed 
surface  t(0)   is  not  parallelizable.  Here  t(0)   is  a 
2-diraensional  torus  having  1-dimensional  invariant  sub- 
manifolds  which  intersect  at  stationary  points  (Fig.  13). 

Gi  =  -ki  sin  01  +  iiFi(y,ei,e2,n) 
02  =  "k2  sin  ©2  +  ^F2(y,9l  ,82,1a) 
y  =  -y  -^   piGCy,ej.,02,|i) 

where  k.  >  0  are  constants,   y  is  distance  measured  along 
the  normal  to  the  surface  and  F.   and  G  e  C°°  in  all 
arguments. 

In  this  example  .3=1  and  X  =  -  maxCki,k2). 
Applying  Theorem  8  we  see  that  jj"  -\  <  l/r ,  r  an  integer 
>  2,   then  for  \j.     sufficiently  small  there  exists  an 
invariant  manifold  y  ='  ^i^i  ,Q2,\a)     where  ^     has  r-1 
Lipschitz  continuous  derivatives  with  respect  to  9;^  ^^^  ©2  and 
is  continuous  in  ^  with  tjiC©!  ,©2  )0)  =  0. 


^ 
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